arXiv:1508.00569v3 [hep-ph] 21 Dec 2015 


Direct and indirect signals of natural 
composite Higgs models 

Christoph Niehoff, Peter Stangl, and David M. Straub* 

Excellence Cluster Universe, TUM, Boltzmannstr. 2, 85748 Garching, Germany 


We present a comprehensive numerical analysis of a four-dimensional model 
with the Higgs as a composite pseudo-Nambu-Goldstone boson that features 
a calculable Higgs potential and protective custodial and flavour symme¬ 
tries to reduce electroweak fine-tuning. We employ a novel numerical tech¬ 
nique that allows us for the first time to study constraints from radiative 
electroweak symmetry breaking, Higgs physics, electroweak precision tests, 
flavour physics, and direct LHC bounds on fermion and vector boson reso¬ 
nances in a single framework. We consider four different flavour symmetries 
in the composite sector, one of which we show to not be viable anymore in 
view of strong precision constraints. In the other cases, all constraints can 
be passed with a sub-percent electroweak fine-tuning. The models can ex¬ 
plain the excesses recently observed in WW, WZ,Wh and resonance 
searches by ATLAS and CMS and the anomalies in angular observables and 
branching ratios of rare semi-leptonic B decays observed by LHCb. Solving 
the B physics anomalies predicts the presence of a dijet or ti resonance around 
1 TeV just below the sensitivity of LHC run 1. We discuss the prospects to 
probe the models at run 2 of the LHC. As a side product, we identify several 
gaps in the searches for vector-like quarks at hadron colliders, that could be 
closed by reanalyzing existing LHC data. 
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1. Introduction 

The Standard Model Higgs boson faces a severe naturalness problem since the presence 
of heavy states associated to a more fundamental theory would lead to enormous correc¬ 
tions to its mass, requiring an extreme fine-tuning to explain the observed value. This 
conundrum can be solved if the Higgs boson is not an elementary scalar but a bound 
state of some new strong interaction. The lightness of the composite Higgs with respect 
to the - as yet unobserved - composite resonances finds a natural explanation if the Higgs 
is a pseudo-Nambu-Goldstone boson (pNGB) of an approximate global symmetry of the 
strong sector [1,2]. To avoid the flavour problems of technicolour theories, the mechanism 
of partial compositeness can be invoked to generate the masses of the SM particles [3]. 


2 
















This mechanism is closely related to the “geometric” generation of fermion mass hierar¬ 
chies from wave function overlaps in models with warped extra dimensions [4-7] , and in 
fact much of the progress in composite Higgs models in the last decades has been made 
using holographic models [8-10]. In these models, the Higgs potential becomes calcula¬ 
ble, leading to additional predictivity. Recently, also four-dimensional models have been 
constructed where the Higgs potential is calculable at one-loop level [11-13]^. These 
models have the advantage that they are simpler in structure than the 5D theories, but 
more general, as they need not necessarily be the low-energy limit of a 5D holographic 
theory. Our aim in this paper is to study one particular implementation of the 4D pNGB 
Higgs, taking into account all relevant experimental constraints. 

Direct and indirect constraints on composite pNGB Higgs models have already been 
discussed in the literature (for recent analyses see e.g. [16,17] for electroweak precision 
tests, [18-20] for flavour physics, [21-25] for Higgs physics, [23, 26-34] for quark part¬ 
ner searches, and [27,35-39] for vector resonance searches). However, a complete and 
simultaneous numerical analysis of all relevant constraints on a particular model is still 
lacking. This is due in part to the fact that a mere parameter scan, as is typically 
done in supersymmetric extensions of the SM, is not feasible since the parameter space 
does not “factorize” into Standard Model (SM) and new physics (NP) parameters; due 
to partial compositeness, the masses of SM particles and the angles and phase of the 
Cabibbo-Kobayashi-Maskawa (GKM) matrix are non-trivial functions of many model 
parameters. As a consequence, finding viable parameter points from a random set of 
model parameters becomes untractable. This problem becomes even more severe once 
the Higgs potential is taken into account, as the Higgs mass and VEV often arise from an 
interplay between gauge and fermion loops which again depends on many parameters. 
For these reasons, numerical analyses of composite Higgs models often have to rely on 
simplifying assumptions, e.g. only considering third generation fermions and their part¬ 
ners - which does not allow including flavour constraints, for instance. Full numerical 
studies of indirect constraints have been performed in warped extra dimensional models 
(without a pNGB Higgs) by making use of approximate analytical expressions for the SM 
parameters [40-45], but this only works for particular representations of the additional 
fermions^. We have overcome these problems by generalising a numerical method first 
proposed in [46] and with the help of a high-performance computing cluster. This allows 
us for the first time to scrutinize one specific model taking into account all relevant 
experimental constraints and to identify novel correlations. 

In selecting a model to analyze in detail, our focus has been to maximize naturalness 
and predictivity, but to be as economic as possible concerning both particle content and 
number of parameters. The model should thus fulfill the following requirements. 

• The symmetry breaking coset should contain custodial symmetry to avoid excessive 
contributions to the T parameter, but no extra Higgs states. This singles out 

^See [14] for an overview of the model landscape and [15] for a comprehensive review of 4D pNGB 
Higgs models. 

^Namely if the left-handed elementary quark doublet mixes with a single composite SU{2)l doublet, 
unlike in the model to be studied below. 
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50(5)/50(4). 

• The Zb]J)L coupling should be custodially protected from tree-level corrections. 
This leaves two possible choices of quark partner representations under the SO{4:) 
symmetry [47]. We choose the one where all quark partners can be embedded in 
two fundamental representations, as in the MCHM5 [48]. 

• The Higgs potential should be calculable. This can be achieved by imposing the 
Weinberg sum rules [13,49]. These are automatically fulfilled in deconstructed 
models like the 4DCHM [11] or the DCHM [12]. We choose the 4DCHM, because 
it features a finite one-loop effective potential already for two sites.^ 

• The contribution to AF = 2 observables, i.e. meson-antimeson mixing, should be 
suppressed compared to the naive anarchic expectation to avoid the ex problem 
[21,41,44]. Several mechanisms have been proposed to address this problem (apart 
from invoking accidental cancellations). We focus on the assumption that the 
composite sector is exactly invariant under a large flavour symmetry which is only 
broken minimally (i.e. by the amount required to reproduce CKM mixing) by the 
composite-elementary mixings.'^ This arguably corresponds to one of the strongest 
assumptions one can make on the flavour structure of partial compositeness, which 
is why we view it as a natural starting point in the search for a model that passes 
precision tests, but is natural in the electroweak sense. 

In summary, we focus on the two-site 4DCHM with quark partners in two fundamen¬ 
tals of SO{5), which we will call M4DCHM5 in the remainder of the paper. For the 
flavour structure of this model, we will consider the four different possibilities studied 
qualitatively already in [18]: an effective U{3)^ [59-61] or U{2)^ [18,61] flavour symme¬ 
try with flavour-invariant composite-elementary mixings either for left- or right-handed 
quarks in both cases, dubbed left- or right-compositeness, respectively. 

We stress that, while we aim to include as many experimental constraints as possible, 
our analysis is on a conceptually different level compared to analyses of weakly-coupled 
renormalizable extensions of the SM, e.g. the MSSM. This is because the models we 
are studying are non-renormalizable with a cutoff in the few-TeV region and contain a 
sector with strong couplings. Consequently, the models are not only less ambitious, but 
also less predictive since contributions from cutoff-scale operators or strong interaction 
effects could potentially spoil the picture obtained from naive computations in the two- 
site picture. Nevertheless, for the observables we are considering, the calculable effects 
often already lead to stringent experimental constraints and we find it unlikely that 
cutoff-scale physics comes to the rescue by cancelling these effects. It should however 
be kept in mind that many of the predictions are afflicted with considerable theoretical 
uncertainties. 


®In the two-site DCHM, a logarithmic divergence spoils the predictivity for the Higgs VEV, but the 
Higgs mass can still be computed [50]. 

^For alternative mechanisms, see [51-58]. 
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2. Model setup 


In this section, we briefly review the M4DCHM and its Lagrangian. For details, the 
reader is referred to the original publication [11], The relation to similar models is 
discussed in [13]. 


2.1. Bosonic part 


The M4DCHM can be understood as a deconstructed description of an extra-dimensional 
Gauge-Higgs-Uniflcation model with a bulk gauge group SO{5) that is broken down by 
boundary conditions on the branes to 50(4) and the SM gauge group. To make the 
model phenomenologically viable, the symmetries are enlarged to include a bulk colour 
sector and an additional U{l)x to match the hypercharge assignments of the SM. So, 
from a 4D point of view, there is a strongly interacting composite sector subject to a 
global symmetry breaking pattern® {SU (3)c x 50(5) x U (l)x )/{SU (3)c x 50(4) x U (l)x) 
and an elementary SM-like sector with gauge group 50(3)^ x SU{2)^ x U{1)^. 

In the two-site model one considers only one level of heavy resonances, thus the spec¬ 
trum contains resonances for the 50(5) as well as heavy gluons and a heavy px^i- 
These resonances mix with their elementary counterparts such that the diagonal group 
becomes the remaining SM gauge group and hypercharge is given as® 

y = T^fl+X (1) 


The bosonic sector of the theory contains the gauge part as well as the sigma model 
describing the global symmetry breaking. 


a 


osonic 


c 


gauge 


+ 


( 2 ) 


The gauge Lagrangian, 

/^gauge = -|tr - ^tr 

- [pG^.uPG^"] - [p^,up^’'] - I PXf.upx^'' 

+ ^ (503 - 9g PGij.y + ^ {90 B^n - gxPXp.)^ , 


(elementary) 
(composite) 
(mixing) (3) 


contains the usual kinetic terms for the elementary SU (3)® x SU (2)® x U (1)® gauge fields, 
where 

IF®^ = IF®;M“^ (4) 

as well as kinetic terms for the gluon-, 50(5)- and f7(l)x-resonances. For the res¬ 
onances it will be useful to group them into SU{2)l, 5[/(2)ij and coset components 
(following [13], the latter we will call “axial resonances” in the following), 

Pm = T“^ + T“« + TT (5) 

®The M4DCHM also contains an additional symmetry breaking (50(5)n x SO{5)r)/SO{5)l+r to 
account for the presence of heavy resonances. 

®See appendix A for our convention for the SO{5) generators. 
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We also introduce explicit mixing terms between the SU{3)- and ?7(l)-resonances with 
their elementary counterparts, which are characterized by the scales fc and fx- 
The sigma model Lagrangian 


C^ = 


fl 


tr 




+ 




55 


( 6 ) 


contains covariant derivatives acting on the sigma model helds bli and ^ 2 , which are 
given as 


= d^VLi - i {go W\ + go'B\ n^ + igp Oi p^, (7) 

VpQ2 = dpQ.2-igpPpi^2- ( 8 ) 


Note that the sigma model fields are uncharged under the global U{l)x symmetry (and, 
of course, they do not carry colour charges). 

We adopt the so-called holographic gauge for the sigma model helds, which is inspired 
by a convenient gauge chosen in the corresponding 5D gauge theory, 


ni(x) =U = exp 


i^(^d{x) T" 

Ji 


^2{x) = l5- 


(9) 


There is also the SM gauge freedom that has to be hxed. Here we adopt the SM 
unitary gauge, such that cra{x) = (0, 0,0,/i(x)). In this gauge the Goldstone matrix 
takes the form 


/I 


U := exp 


i^(Taix)J° 

Jl 


V 


cos 
— sin 


h{x) 

fi 

h(x) 



( 10 ) 


Writing the Lagrangian as above in holographic gauge leads to a mixing term of the 
form 

( 11 ) 


9p fi a4 dph. 

One can get rid of this term by a held redehnition, 

a4 —)• O4 — , 

gp J 2 


i, ^ hr, 

k^jh, 


( 12 ) 


where / is given by := By this transformation the mixing term van¬ 

ishes and the composite Higgs kinetic term is canonically normalized. As a result, all 
dependencies on the Higgs held are given via 


Sh = sin 



(13) 
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2.2. Fermionic part 

For the fermionic part of the model, we distinguish between the quark and the lepton 
part, 

-^fermionic — -^quark “1“ -^lepton- (1^) 

As in the boson sector, the quark Lagrangian contains elementary, composite and mixing 
parts. 


Tquark = iff Lp<f L + iu^ rTPvP R + idP Rpd^ R 

(elementary) 

T comp't^^comp T ^Tcompcomp 

(composite) 



- {mYjj + Yjj) SulSuR - fnYuQuiQuR + h.c. 



+ Au^ {QuR + Sur) + (^QuL + SuL^ + h.c. (mixing) 

+ (u o d) (15) 

Here we have two bidoublets, Q and Q, and two singlets, S and S, for every flavour. 
For the kinetic terms we used an SO{5) notation where we combined the singlets and 
bidoublets into 50(5) fundamentals: 'I'comp = {Q,S)u,(R ^comp = {Q,S)u,d- For these 
the covariant derivatives are then defined as 


F^/^'kcomp — {dfM iQG PGfi ipp iQX 9X PXp)comp (16) 


and the same for 'kcomp- The U{l)x charges are assigned to match the hypercharge of 
the SM. Thus, the fundamentals Tcomp and Tc^mp have = | and = “5- 
The elementary fields are embedded into (incomplete) 50(5) fundamentals via 


^uL = 


1 

71 


^dL — 


1 

71 


/ d, \ / 0 \ 


-idi 


0 

UL 

II 

0 

iUL 


0 

\ 0 y 


\UR J 

( UL \ 


( ^ \ 

iUL 


0 

-dL 

II 

0 

Ul 


0 

\ 0 y 


\ dR / 


(17) 


(18) 


Since we are mainly interested in the interplay between quark flavour measurements 
and the Higgs potential, we do not consider effects of partial lepton compositeness in 
this work. Indeed, if the compositeness of the left- and right-handed lepton chiralities 
are comparable, they are required to be small due to the leptons’ lightness and their 
impact on the observables to be considered below is expected to be small. Moreover, 
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flavour-changing interactions are strongly constrained by negative searches for charged 
lepton flavour violating processes. In practice, we simply consider elementary leptons 
with direct bilinear couplings to the Higgs field, 

/:iepton = ilL^L + ^ (19) 

where, just as for the elementary quarks, the covariant derivatives are understood as 
couplings to the elementary SU{3)^ x SU{2)^ x C/(l)° gauge fields. We note however 
that a significant degree of compositeness for some of the leptons could be motivated 
experimentally, e.g. to reconcile radiative electroweak symmetry breaking (EWSB) with 
naturalness in the absence of light top partners [62] or to explain the hints for violation 
of lepton flavour non-universality in B decays [63] . These effects are beyond the scope 
of our present analysis. 

Let us note here that many of the above model parameters are correlated if they 
originate from an extradimensional gauge theory, e.g. coupling constants are generated 
by overlap integrals of Kaluza-Klein mode functions. In our numerical analysis, we will 
not impose such relations but instead try to be as general as possible to explore the 
viability of purely 4D pNGB Higgs models, regardless of whether a dual 5D description 
exists. 

2.3. Flavour structure 

As noted in the introduction, we assume the strong sector to be invariant under a flavour 
symmetry, only to be broken by the composite-elementary mixings. We consider four 
possibilities (see [18] for a thorough comparison), 

• In U (3)lc (LC for left-handed compositeness), the strong sector is invariant under a 
t/(3) symmetry^ that is broken by the right-handed composite-elementary mixings. 

• In U{3)^q (RC for right-handed compositeness), the strong sector is invariant 
under a 17(3) x 17(3) which is broken by the left-handed composite-elementary 
mixings. 

• The U{2)\^q and U{2)\q models are analogous, but restricted to a smaller symme¬ 
try only acting on the first two generations of composite quarks. 

In the U{3)^ models, the matrices in the composite part of the fermion Lagrangian (the 
3rd and 4th lines in (15)) are proportional to the identity, while in the U{2)^ models, they 
are of the form diag(a, a, b). The main difference is the form of the composite-elementary 
mixings. We give their explicit forms in appendix C. 


^The cube in U{3)^ refers to the fact that after the breaking, the SM quark sector is approximately 
invariant under a U{3)q x U{3)u x U{3)d- 
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2.4. Higgs potential 

The explicit breaking of the global symmetries by the mixings of the composite reso¬ 
nances with the elementary sector generates an effective potential for the NGB field (to 
be identified with the SM Higgs) at the quantum level, so that it acquires a mass and a 
VEV, breaking electroweak symmetry. In general, the effective potential is UV-sensitive 
and not necessarily calculable. In the M4DCHM it is finite at one loop, making the 
model predictive under the assumption that higher loop contributions are subleading 
with respect to the (calculable) one-loop contribution. We will rely on this assumption 
in the following. 

At one loop, the effective potential is given in terms of all the n-point correlation 
functions of the Higgs and therefore contains a gauge as well as a fermion contribution. 
It can be calculated by the Coleman-Weinberg formula [64] 


Ves{h) = (2tr [Mf{h)] h? - tr log [A^ 

+ tr log [M‘^{h)] 


( 20 ) 


where Mf{h) denote the Higgs-dependent mass-(mixing)-matrices^ which we give in 
appendix B and 

{ 3 for neutral gauge bosons 
6 for charged gauge bosons 
— 12 for (coloured) Dirac fermions 

Here we explicitly showed the dependence on the cutoff of the theory. For a non- 
renormalizable effective theory these UV dependent terms spoil the predictivity, such 
that one has to demand the following relations for ensuring the calculability of the Higgs 
potential: 


tT[Mf{h)] -tv[Mf{h = 0)]=0, (21) 

tr [iMhh)f] - tr [{M^h = 0))^] = 0, (22) 


where it was taken into account that the constant term of the potential is not physical. 
These relations are just a reformulation of the Weinberg sum rules of the fermion and 
gauge sector that are usually imposed to guarantee a finite potential [13,49]. For the 
quark sector, these relations represent a generalisation of the Weinberg sum rules to the 
three family case. 

In deconstructed models the Higgs potential is usually protected by the higher dimen¬ 
sional gauge symmetry, such that the Weinberg sum rules are automatically satisfied. 
This is the case for the M4DCHM as well [11,13]. Note that in this case also the scale 
dependence cancels from the effective potential. Then the expression for the potential 
simplifies to 

Vesih) = Y ^ogim^ih)), (23) 

all particles 

®For fermions these are given by M^{h) = M{h)^M{h). 
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Figure 1: Notation for Feynman rules used in this work. 


where mi{h) denote the masses in the mass basis, i.e. the singular values of the mass 
matrices. 


3. Experimental constraints 


In this section, we discuss all the experimental constraints that we impose in our anal¬ 
ysis. Since approximate analytical expressions for most of the observables have already 
been provided elsewhere (see in particular [18,19,46,57]), we focus on discussing the 
numerical computation and on specifying our treatment of theoretical and experimental 
uncertainties. Section 3.1 specifies how we compute the masses and couplings of the SM 
states, including among others the Higgs mass and VEV as well as the CKM matrix, sec¬ 
tion 3.2 discusses indirect constraints, including electroweak precision observables and 
flavour physics, while section 3.3 deals with the direct bounds on fermion and vector 
resonances. 

Since the numerical computation of the observables involves masses and couplings in 
the mass eigenstate basis (mb), we fix our notation for the couplings by specifying the 
Feynman rules in figure 1. In the gauge basis (gb) the Lagrangian contains non-diagonal 
mass matrices for vector bosons and fermions as well as interaction terms connecting 
both kinds of fields, which are (schematically) given as follows 


After EWSB one can go to the mass basis by unitary transformations. 


4 ' 


(gb) 
L,R i 




iL,R) 






j 


(mb) 
L,R i’ 


A 


(gb) 

/i i 


= [b,|y A^, 


(mb) 


such that the couplings as defined in figure 1 can be calculated as 


(mb) 


= IbJla. 


- abc 


dgb) 


J ijk 


(24) 

(25) 


(26) 
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3.1. Standard model masses and couplings 

3.1.1. Higgs VEV and masses of SM states 

The tree-level masses for fermions and gauge bosons are obtained by diagonalizing the 
mass matrices given in appendix B after EWSB. Since the interaction terms of the 
Goldstone bosons generate mixing terms between elementary fields and the composite 
resonances, all masses of fields with SM-like quantum numbers will depend on the VEV 
taken by the pNGB, i.e. they depend on = sm((h} jf). In the model used in this 
work this quantity is not a free parameter, but it is calculable as the minimum of the 
loop-generated effective potential. 

In practice, we calculate the s/^-dependent masses of all particles and use (23) to 
calculate the effective Higgs potential. Then, the correct value of is obtained by nu¬ 
merically minimizing the potential. We also explicitly demand that the found minimum 
is non-trivial, otherwise we discard the parameter point. 

The value of is fixed in our numerical analysis by imposing the tree-level value of 
the Fermi constant in muon decay as a constraint, 

^tree _ ^ ^ _ 

" “ “ V2(siy-‘p • 

which is valid up to negligible vector resonance exchange contributions. Since we only in¬ 
clude the tree-level we add a relative theoretical uncertainty of 1%. Loop corrections 

to Gf^ are effectively included in terms of the T parameter, see section 3.2.1. 

Once is known, the Higgs mass can be calculated as the curvature of the effective 
potential at its minimum, 

= ^2 ^ ■ (28) 

For the W, Z, and top masses, we directly interpret the masses obtained from diag¬ 
onalizing the mass matrices as MS running masses at the scale mt- We add a relative 
theory uncertainty of 5% to account for this crude assumption. In principle, we could 
compute the one-loop matching corrections to the masses to get a more reliable estimate. 
In practice, this is not feasible because the composite-elementary mixing means that the 
numerical computation of a large self-energy matrix, e.g. 27 x 27 in the case of the top 
quark, would be necessary which quickly leads to excessive computing times. 

For the light quark masses, we also interpret the tree-level masses as MS running 
masses at mt] then we use RunDec [65] to run them to the relevant scales where they 
can be compared to the PDG averages [66]. 

3.1.2. CKM matrix 

In the SM, GKM elements are determined from a global fit to weak decays mediated by 
tree-level W exchange as well as loop-induced meson-antimeson mixing observables. In 
the presence of NP, the latter are susceptible to NP contributions as will be discussed 


(27) 
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below in section 3.2.4. But even the tree-level processes receive corrections in a composite 
Higgs framework that lead to relevant constraints. The reason is that the 3x3 quark 
mixing matrix is no longer unitary in the presence of composite-elementary mixing, but 
becomes part of a larger (27 x 27) mixing matrix among quarks and quark partners. 
Deviations from CKM unitarity, predicted by the SM, can thus be used to constrain 
quark compositeness. 

To compare to the absolute values of CKM elements measured in experiments, one 
can define effective CKM elements from ratios of W couplings, 

IdWuidj 

The \Vij\ obtained in this way can be directly compared to the elements extracted in 
experiments assuming the SM as long as right-handed W couplings and contributions 
from tree-level heavy resonance exchange can be neglected. We do take these two effects 
into account in our numerics, although they turn out to be negligible. 

In our numerical analysis, we include five CKM elements that are directly measured 
in tree-level processes. 


(29) 


\Vud\ from superallowed nuclear beta decays, 

|I4s| from K —)• -kIv decays, 

\Vuh\ from inclusive B —)• and exclusive B —)• -irii' decays, 

\Vcb\ from inclusive B —?• Xcii' and exclusive B —)• D*iiy decays, 

\Vtb\ from the cross-section of f-channel single top production at LHC. 


The measured values and references are given in table 1. In the case of |I4b| and |I4fe|, 
there are long standing discrepancies between the determinations from inclusive vs. ex¬ 
clusive B decays. Since these tensions cannot be resolved in our model, we use the PDG 
prescription [ 66 ] to rescale the discrepant measurements. We multiply the uncertainties 
of \Vub\ by a factor of 1.9 compared to the ones given in table 1, and a factor of 2.9 in 
the case of \Vcb\- 

\Vud\ and |I4s| are important because in the SM, they are constrained by the unitarity 
condition on the first row of the CKM matrix. 


1 = \Vud? + IK.P + \Vub\^ ~ iKdl" + IK.P , (30) 


where \Vub\ is numerically negligible. The smallness of \ Vub\ and |I4fc| is also why, in the 
SM, \Vtb\ ~ 1 holds up to a permille level correction. Partial compositeness can lead to 
a deviation from both relations (see e.g. [16,18,60]). 

Finally, we also include the CKM angle 7 that is measured via the interference of 
b —>■ cus and b —>■ ucs amplitudes in H —)> DK decays. Again in the case where right- 
handed W couplings and direct vector resonance contributions can be neglected, 7 can 
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Rb 

0.21629(66) 

[67] 

\Vud\ 

0.97417(21) 

[68] 

Rc 

0.1721(30) 

[67] 

\Vus\ 

0.2249(8) 

[69] 

Rh 

20.804(50) 

[67] 

1 

(3.72 ±0.16) X 10-3 

[70] 

QQ 

9ww 

0.86 ±0.17 

[71,72] 

1 I^blin 

(4.33 ±0.28) X 10-3 

[73] 

QQ 

9zz 

1.18 ±0.39 

[71,72] 

1 ^cb 1 ex 

(3.904 ±0.075) X 10-^ 

[74] 

QQ 

9gg 

1.12 ±0.22 

[71,72] 

\Vcb\in 

(4.221 ±0.078) X 10-2 

[75] 

gg 

r-T+T- 

0.97 ±0.39 

[71,72] 

\Vtb\ 

0.998 ±0.041 

[76] 

AMk 

3.483(6) X 10-15 GeV 

[66] 

\^k\ 

2.228(11) X 10-3 

[66] 

AMd 

0.510(3)ps-i 

[77] 

^IpKg 

0.682(19) 

[77] 

AMs 

17.761(22) ps-i 

[77] 

(ps 

-0.010 ±0.039 

[78] 

S 

0.05 ±0.11 

[79] 

7 

(72.9 ±6.7)° 

[80] 

T 

0.09 ±0.13 

[79] 





Table 1: Values of the experimental constraints used in the numerical analysis. For details and 
the treatment of theoretical uncertainties, see main text. 


be computed from the tree-level W couplings as 

ffWud 9wub ^ 
9wcd 9wcb / 


7 = arg 


(31) 


This expression is independent of phase conventions. For the experimental value in 
tab. 1, we symmetrize the value obtained by the CKMfitter collaboration from a fit to 
all experiments. 


3.2. Indirect constraints 
3.2.1. S and T parameters 

By construction, the T parameter does not receive a contribution at tree level in pNGB 
models based on the 50(5)/50(4) coset. At one loop, the dominant contribution typ¬ 
ically comes from fermion loops involving, in particular, the top partners. In addition, 
the modification of the gauge boson couplings to the Higgs and the electroweak would- 
be Goldstone bosons leads to an “infrared-log” contribution [81,82]. Finally, also loops 
involving the heavy spin-1 resonances can contribute (see [17] for a recent discussion). 
For simplicity, in our analysis we restrict ourselves to the fermion contribution, which is 
finite and gauge-independent. It can be computed numerically as 




TT^ 

rr _ ‘-'-WW 
~ ^2 
W 


m 


niz 

2 

m% 


(32) 


where the masses are tree-level masses, and 
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fijj 


2 2 


dVfifj? + \9vfifj?'') 

+ ^^[gvufjdvuf) 


( 33 ) 


is the fermion contribution to the transverse part of the vacuum polarization. The sum 
runs over all SM fermions and quark resonances. The Passarino-Veltman function is 
defined as in [83] and the function H can be found e.g. in [84]. 

In contrast to T, the S parameter arises already at tree level, effectively leading to a 
lower bound on the mass scale of the spin-1 resonances. In models where T = TgM at 
tree level, the NP contribution to S can be obtained numerically as 


Clem'S* It=0 — ^ i.^W ^eff) i 


(34) 


where = 1 — and the effective weak mixing angle is defined via the leptonic 

forward-backward asymmetry, 


^ 9zee + ffZee 
9zee ~ 9Zee 


sin^ ^eff 


1 + X 
4 


(35) 


Experimentally, a recent global fit of electroweak precision data finds [79] 


5 = 0.05 ±0.11, 


T = 0.09 ±0.13, (36) 


with a correlation coefficient of ±0.9. Since we neglect gauge contributions to T and 
all loop contributions to S, in our numerical analysis we further assume uncorrelated 
theory uncertainties of 0.05 for S and 0.10 for T, which we combine with the correlated 
experimental uncertainties. The size of these theory uncertainties is chosen to encompass 
the typical size of the neglected “IR-log” contributions to S and T. 


3.2.2. Z decays 

Due to the large degree of compositeness required for the left-handed top quark (and 
thus also b quark), the partial width of the Z into b quarks measured at LEP provides 
a powerful constraint on models with partial compositeness. While our model features 
a custodial protection of this coupling, the observable is still important to constrain the 
subleading composite-elementary mixing of the b^. In the flavour-symmetric models, 
also the partial widths into lighter quarks lead to constraints. We include the following 
observables in our analysis. 


_ nz ^ bb) _ T{Z ^ cc) _ TjZ ^ qq) 

T{Z ^qq)' T{Z ^ qq)' ^ T{Z ft) 


(37) 


where r(Z —)• qq) implies a sum over all quarks but the top. We compute only the 
tree-level corrections at zero momentum to these observables (see [16] for a discussion of 
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effects beyond this limit). We add the higher-order SM contributions (see [85]) numeri¬ 
cally to reproduce the correct SM predictions in the absence of NP contributions. The 
experimental measurements are listed in table 1. 

A comment is in order on the loop corrections to Z —)• 6^6^, which we have not 
taken into account. Although corrections are already generated at tree level in the 
M4DCHM5, these are suppressed by the custodial protection mechanism, which however 
is not active at loop level. In [81,86,87] it was shown that in similar models as the ones 
we are studying, there is a correlation between fermionic loop corrections to the T 
parameter and the loop correction to Z —?■ 6^6^. For a heavy new physics scale, this can 
be understood as being due to renormalization-group mixing of dimension-6 operators 
invariant under the SM gauge symmetries. Considering the operators (in the notation 
of [88] and in the basis where the down-type quark mass matrix is diagonal) 

( 937 ^ 93 ), Q^ll = (<(>^*^J)0)(g3'rV93), (38) 

the correction to the left-handed Z coupling to bottom quarks satisfies 

^9zhh ot ■ (39) 

Custodial protection implies ^ subleading mixing effects, implying a 

vanishing correction at the matching scale. Since the two Wilson coefficients run 

differently [89-91], a non-zero correction is induced at the electroweak scale which is 
proportional to the matching scale value of The quantum corrections leading to 

this running induce at the same time a non-zero T parameter which is also proportional to 
the matching scale value of and thus correlated to Sgzbb- ^ positive contribution 
to the T parameter, the sign of this correlation leads to a negative contribution to Rb 
that is disfavoured by experiment [81,86,87]. Thus, parameter points with a large 
positive contribution to the T parameter might be excluded by taking into account the 
one-loop corrections to ^gzbb- ^ challenge of taking this loop contribution into account 
is that it involves Passarino-Veltman functions at non-zero external momentum with 
three propagators. Due to the large number of states in the M4DCHM5, this would 
significantly increase computing time, so we are not able to take this effect into account. 
It should thus be kept in mind that our results might be optimistic in the sense that 
we might keep points that are possibly excluded. A dedicated analysis of the impact of 
higher order corrections to Sg^^^ would be worthwhile. We also note that the tension 
between the constraints on the T parameter and Sg^bb might be relaxed by including an 
additional level of resonances [87]. 

3.2.3. Higgs production and decay 

We compute the modification of the Higgs partial widths rx = r(/i —)■ X)/T{h —)• Ar)sM 
at tree level for X = WW, ZZ, bb, and and at one-loop level for X = gg and 

® The Wilson coefficient C^u associated with the Operator D ^(^)(m 37 ^M 3 ) that also enters 

the RGE induced contributions to vanishes due to custodial protection. 
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77 - We take into account the loop contributions from all SM and heavy fermions and 
vector bosons. The signal strength in a particular final state, assuming pure gluon fusion 
production, can then be obtained as 


= 


rxrgg 

^tot 


(40) 


where rtot = is the modification of the total width. 

We use ATLAS and CMS measurements to constrain the signal strength. In the case 
of ATLAS, the gluon fusion result is given explicitly. In the case of CMS, we use the 
“0/1 jet” result for WW, ZZ, and r+r”, and the “untagged” result for 77 . We naively 
combine the ATLAS and CMS results for each final state, using the PDG prescription 
to enlarge the error in the case of poor agreement. The resulting constraints are listed in 
table 1. We neglect the correlations between individual measurements. Since the /i —)• 66 
signal strength is only measured in the case of vector boson associated production, we 
do not include it in our numerical analysis. 


3.2.4. Meson-antimeson mixing 

The meson-antimeson mixing amplitude for the neutral meson (= K^, Bg, or 
D^) can be written as 


with q = u or d. The loop-induced SM contribution is discussed e.g. in [92]. The sum 
contains NP contributions due to tree-level vector resonance exchange that are encoded 
in the Wilson coefficients of the following AT = 2 operators. 


Qv'lR ~ {QLl^QL)idRl^QR) ) Q'sLR — iQRQL)iQLQji) 5 

that can be written in terms of the Feynman rules defined in figure 1 as 


\ ^ 

'-^VLL — r, / , 




/ L \ 2 
1 / dpoqjqk \ 


m 


PG 


(42) 

(43) 


(44) 


/^QiQj _ ^QiQj I 

^VRR - ^VLLIl^R 


^qiqj __ \ ^ 

'^VLR — 


g\ 9^0 

^p’qjQk^Piqjqk 

2 

m 0 
Pi 


I 1 dpaqjqkdpoqjqk 

' ^ o 


6 


m 


PG 


(45) 

(46) 


L/ R 

f^qiQj _ ^PGqjPkdpGqjqk 
'^SLR ~ ^2 

"'PG 


(47) 


These expressions are valid at the matching scale of NP and SM, while (41) depends on 
the values at the hadronic scale /i; that is chosen conventionally as mb for B^^s mixing. 
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3 GeV for mixing, and 2 GeV for mixing. In principle, the correct matching 
scale is set by the mass scale of the heavy resonances. However, in our numerical scan, 
we often encounter vastly different scales for the different resonances. Gonsequently, we 
have decided to simply match the tree-level Wilson coefficients to the SM at the scale 
mt and to neglect the RG evolution above m*. A more complete treatment, including 
intermediate thresholds, is beyond the scope of our present work. We note that the 
RG evolution typically makes the NP effects larger at low scales. In that sense, our 
treatment leads to more conservative bounds. All relevant anomalous dimensions for 
the evolution below mt can be found in [93]. 

The matrix elements in (41) depend on meson decay constants and bag parameters, 
both of which can be determined from lattice QCD. They can be written as 

, (48) 


where 


Bvll — Bvrr “ 3-^1 ( w ) 


nM _ ^ 

‘^VLR ~ I I 

6 \mg^+mg^ 




nM _ 
^SLR — 7 


rriM 


mq- + m, 




<(w). 


(49) 

(50) 


For the lattice predictions of the decay constants as well as the bag parameters Bi for 
Bfi and Bg mixing, we use ref. [69], for the kaon bag parameters we use ref. [94], and for 
the charm bag parameters ref. [95]. 

We use the following observables sensitive to NP in the meson-antimeson mixing am¬ 
plitude. 


• The mass differences in the B^ and Bg systems, 

AMd,g = 


(51) 


For the theoretical uncertainties, which are dominated by lattice uncertainties, we 
take 10.2% relative uncertainty for AM^ and 7.6% for AMg. Note that we do not 
have to account for uncertainties due to GKM elements as these are allowed to 
vary in our scan. We further take these lattice uncertainties to be correlated with 
a coefficient of 0.17, since the ratio of the relevant lattice parameters is known 
more precisely than for the individual systems. 

• The mixing-induced CP asymmetry in B^ —)• J/^Kg, 

Si.Ks = sin (avg(^M^2‘')) > (52) 


which in the SM measures sin 2/3. We add a theory uncertainty of 0.01 to account 
for possible penguin pollution [96]. 
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• The sine of the Bg mixing phase as obtained from an average of the mixing-induced 
CP asymmetries in Bg —>■ J/il^ K~^K~ and Bg —)• decays^®, 

sin(/)s = sin ^arg^Mj^®^^ . (53) 


In this case, we add a theory uncertainty due to penguin pollution of 0.017 [96]. 
The parameter for indirect CP violation in mixing, 


Im M 


ex = Ke 




y/2A.Mx 


(54) 


where the experimental value for AMk can be used. For the (non-CKM) theory 
uncertainty on |ex|, we take a relative uncertainty of 11%. 


The mass difference in mixing. 


AMk = 2Re( 


(55) 


The SM contribution to AMk is plagued by large uncertainties due to long-distance 
contributions. Although first results are available from lattice calculations [99], 
these are still for unphysical kinematics. Thus we conservatively allow the NP 
contribution to saturate the experimental central value at la (i.e. at 3a, we allow 
points where the NP contribution is three times the experimental central value, 
implying a necessary cancellation with the SM contribution). 

We do not impose mixing observables as constraints, as they are expected to receive 
small NP contributions in models with minimally broken U{2)^ [61], but we will discuss 
predictions for them in section 4 below. 


3.2.5. Rare B decays 


The b ^ s'j transition arises first at the one-loop level; approximate analytical results 
as well as generic formulae that can be used in a numerical analysis have been presented 
in [19]. We include the constraint from the B —)• Xgj branching ratio, that agrees well 
between SM prediction [100] and experimental world average [77] , 

BR(R ^ Xg'y)sM = (3.36 ± 0.23) x 10"^ , (56) 

BR(R ^ Xg'y)wA = (3.43 ± 0.22) x 10“^ . (57) 

For the NP contribution, we use the following formula (cf. [101]), 


BR{B -A Xg-f) 
BR{B -A Ws7)gM 


1 

(|Cf(m6)SM|2 + iV^) 


(|C'7®(?7lfe)P + |C'7(w-b)P + ATy)^ , 


(58) 


^°Here we have used the average performed by the LHCb collaboration. Very recently, a measurement 
with comparable precision has been presented by CMS [97]. The observable has also been measured 
by ATLAS [98]. 
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where = 3.6 x 10“^. 

The imaginary part of the Wilson coefficients and the relative size of the left- and 
right-handed Wilson coefficients can be constrained by other processes, most notably 
B —> K *angular observables. We do not impose these additional observables as 
constraints, but will discuss predictions for them in section 4. 

NP contributions to semi-leptonic FCNC decays of B and K mesons arise already at 
tree level, mediated by the Z boson that can obtain flavour-changing couplings to quarks 
at tree level as well as by heavy neutral vector resonances (for a thorough discussion of 
these effects in composite Higgs models, see [46]. Similar effects are obtained in models 
with a warped extra dimension [42,44]). Writing the four-fermion operators as 

Qvab = , (59) 


with A,B = L, R, the Wilson coefficients are obtained in analogy with section 3.2.4 as 


„didjU _ 

^VAB — 


9zdjdk9zii 


rUr 




(60) 


Here we have explicitly included the Z contribution as it contributes formally at the same 
level in vjf as the heavy resonance exchanges. The smallness of the flavour-changing 
coupling (which only arises after EWSB and is of order n^//^) is compensated by the 
absence of the suppression by the heavy resonance mass in the propagator. 

One can map the coefficients (60) onto the traditional operator basis for di -A- 
transitions as 


^9 — ^Hj y^VLR + ^VLL J 

dj^dj 1 2 f ^d{dj££ ^d{dj^£\ 

^9 — y^VRR ^VRL ) 


*-^10 ~ y^VLR ~ ^VLL J 

^/didj _ A 2 /^didjU _ ^didjU\ 

*-^10 “ y-'VRR ^VRL ) 


(61) 

(62) 


where 

A?. = —- - -. (63) 

V2GFaemVtiV* 

The primed coefficients are only generated at a very suppressed level in the flavour- 
symmetric models we consider. Since we are assuming leptons to be elementary, all 
Wilson coefficients are lepton flavour universal. Relaxing this assumption, the recent 
hint for lepton flavour non-universality can potentially be explained as well [63], but we 
are not considering this possibility here. Since the lepton-Z couplings are SM-like to 
an excellent precision in our setup, the Z-mediated contributions fulfill the well-known 
relation Cg = (4s^ — l)C'io, i-e. they mostly contribute to Ciq. 

Concerning the resonance-mediated contributions, as mentioned above, they formally 
contribute at the same order in u// as the Z contributions. Their couplings to elementary 
leptons however only arise through mixing of the composite and elementary vectors, 
so the resonance-mediated contributions are expected to be parametrically suppressed 
compared to the Z-mediated ones by a factor g\l where 5ei,co are generic elementary 
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and composite gauge couplings. Still, there are parts of parameter space where these 
contributions can be relevant. To understand their impact, it is instructive to work in 
a basis where instead of the three electrically neutral electroweak resonances p^, p^, 
and p^, one works with three linear combinations that, before EWSB, couple to the 
same quantum numbers as the Z, the photon, and one which does not couple to the 
leptons at all. The first two states are analogous to the KK Z and the KK photon 
in Randall-Sundrum models (cf. [43,102]). This basis is relevant because the “KK Z” 
contribution leads to Cg = (4s^ — l)C'io just as the Z-mediated one, while the “KK 
photon” contribution affects only Cg and not Cig. In addition, the part of the “KK Z” 
contribution that involves the composite-elementary mixings is forbidden by the 
same custodial protection that protects the ZhibL coupling, while a similar protection 
is absent for the “KK photon” contribution. This is particularly relevant in C/(3)^(-., 
where the correction involving is flavour diagonal in the mass basis [18], cf. (90). 
As a consequence, the Z-mediated as well as the “KK Z” contribution to the AT = 1 
operators vanish, while the “KK photon” contribution can be nonzero. 

Recently, a number of tensions between measurements and SM expectations have 
appeared in several observables in rare 6 —)• s decays. This includes in particular 

• A suppression of the angular observable P 5 in R —)• K* [103-105]; 

• A suppression of the branching ratio of Bg —)• (pp'^p~ [105,106]; 

• A suppression of Rk, the ratio of branching ratios of B~^ —)• K~^p'^p~ and B~^ —?■ 
K+e+e- [107]. 

While the first two of these anomalies could be due to unexpected hadronic effects 
(see e.g. [108,109]) and the last one due to a statistical fluctuation, all of them could 
be explained consistently by a negative NP contribution to the Wilson coefficient Cg® 
(a positive contribution to CJ’q is allowed in addition) with muons only [110-115]. In 
composite Higgs models, such lepton flavour non-universal contribution was shown by 
us to arise if muons carry a significant degree of compositeness [63]^^. In the present 
setup, since we are considering elementary leptons only, all effects are lepton flavour 
universal. Although in this case, the deviation in Rk cannot be explained, the overall 
agreement with the data could still be significantly improved compared to the SM if 
there are (lepton flavour universal) NP contributions in Cg® (and possibly C]’q), which 
can resolve the tensions in R —)> K*p~^p~ angular observables and various branching 
ratios and give a good fit to the data [115]. As discussed above, such contribution can 
arise from “KK photon”-like resonance exchanges. 

In view of these tensions, we do not impose semi-leptonic 6 —>• s transitions as con¬ 
straints in our numerical analysis, but rather consider the predictions for them a poste¬ 
riori. 

We do however include the branching ratio of Bg —)• p'^p~ as a constraint. This 
branching ratio, which has reduced theoretical uncertainties compared to the semi- 
leptonic decays, was recently observed by LHCb and CMS [117] in agreement with the 

Another possibility is to introduce composite leptoquarks [116]. 
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SM expectation [118], 


BR{Bs //+/i“)sM = (3.65 ± 0.23) x 10“® , (64) 

BR{Bs ^ //+/i-)exp = (2.8^°;^) x IQ-^ . (65) 

In the presence of new physics, the branching ratio is modified as 

BR(R, ^/i+/i-)sM |(C?g)sMP ■ ^ ’ 

Again, the imaginary parts and the chirality structure can be constrained by other 
observables in exclusive and inclusive semi-leptonic decays. 


3.2.6. Contact interactions 


Four-quark contact interactions are constrained by measurements of the dijet angular 
distribution at LHC. These constraints become relevant when some of the first-generation 
quark fields have a signihcant degree of compositeness. This is unavoidable in the U{3)^ 
models, but also occurs in part of the parameter space of the U{2)^ models. The relevant 
four-quark operators involve only the first generation quarks as the contribution from the 
other generations is PDF-suppressed. The Wilson coefficients are computed analogously 
to the AF = 2 ones in section 3.2.4. Experimental bounds are usually quoted on 
operators in an SU{2) x U{1)y gauge-invariant basis. Using the notation of [119], their 
Wilson coefficients can be related to the ones in the low-energy basis as 


cS = C\ 

eg = c; 


uu 

VLL 


uu 

VLR 


.(1) 


_ /-luudd 

— Lxi 


~'uudd ^VRR 
^( 1 ) _ r^uu 


1 ^ /-luddu 

+ q^VLL 5 

^qq 

/^UU 

— ^VLL ; 

( 67 ) 

^ /^UU 

~ q^SLR 5 

^qu 

^ y^UU 

— -q^SLR 5 

( 68 ) 

1 ^ ^uddu 

+ -^^VRR 5 

uudd 

ry^uddu 
—^^VRR 5 

( 69 ) 

1 

^dd 

y^dd 

— ^VRR •> 

( 70 ) 


fl 8^ 

and with the appropriate replacement u ^ d for . In addition to the operators in 
(42), (43) and the ones with d —)> n, we have defined 

QvLL = {uLl^dL){dLrUL) , Qr" = {uLrUL){dLl^dL) , (71) 


as well as L —)• R. 

The Wilson coefficients of the four-quark operators in the low-energy basis can be 
computed analogously to the AF = 2 Wilson coefficients in section 3.2.4. However, 
an important difference is that the measurement of the dijet angular distribution at 
LHC involves processes at much higher energies compared to meson decays. The EFT 
approach is only valid if the exchanged resonances are much heavier than the typical 
energy scale of the process in question. In [120], it has been shown that for resonance 
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masses below about 5 TeV, the contact interaction bounds become much weaker than a 
naive application of the EFT would suggest. To account for this fact in an approximate 
way, we follow a prescription advocated in this paper and multiply every individual 
contribution to the four-quark operators arising from exchange of a resonance with mass 
nip. by a correction factor (1 + j, adopting C = 1.3 TeV as a rough estimate 
based on a numerical analysis of the full mass dependence in two benchmark scenarios 
[ 120 ]. 

ATLAS and CMS have presented constraints on contact interactions using the full run- 
1 data set. However, the constraints are only quoted for a single operator (in the case of 
ATLAS) or for individual operators, but only allowing one at a time (in the case of CMS). 
In our case, multiple operators might be present simultaneously, and the operators with 
right-handed quarks typically differ for up- and down-type quarks. The full dependence 
of the dijet angular distribution on all operators has been discussed in [119] and simple 
formulae for the impact of the operators in specific rapidity bins have been presented 
there for the 7 TeV LHC. We use these results, updated to the 8 TeV LHC, to obtain 
the relative contributions of individual operators to the differential cross section, while 
we use the bound on the Wilson coefficient c^qq quoted by the experimentalists for the 
normalization. Details on the procedure are discussed in appendix D. 

3.3. Direct searches 

In addition to the indirect searches, i.e. precision constraints from flavour, electroweak, 
and Higgs physics, composite Higgs models are also subject to increasingly strong direct 
constraints from searches for composite resonances at the LHC. Since we are focusing 
on a model with a minimal Higgs sector and we are ignoring partial compositeness of 
leptons, in our case the relevant searches are the ones for quark partners, to be discussed 
in section 3.3.1, and for spin-1 resonances, to be discussed in section 3.3.2. 

3.3.1. Quark partners 

Pair production of heavy quarks and subsequent decay to SM quarks and weak bosons 
has been searched for at Tevatron and LHC. Recently, also final states involving the 
Higgs have been included in the searches. In the simplest case where only decays to third 
generation quarks and a W, Z or Higgs are allowed, these channels can be combined to 
obtain stringent bounds on the masses of vector-like 3rd generation quark partners. In 
our numerical analysis, we aim to be more general since in principle, a quark partner can 
have several relevant decay modes involving SM or partner quarks, 3rd or light generation 
quarks. To this end, we compute the production cross section times branching ratio of 
each quark partner in each experimentally relevant decay mode, and compare it directly 
to the upper bounds on this quantity provided in the experimental analysis. 

For the pair production cross section, we simply take the model-independent NNLO 
QCD production cross section for a heavy quark computed in Hathor [141]. This means 
we neglect 
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Decay 

Experiment 

[TeV] 

Luminosity [fb 

Analysis 


Q^tW 

CMS 

7 

5 

B2G-12-004 

[121] 


ATLAS 

7 

1.04 

EXOT-2011-28 

[122] 


CDF 

1.96 

4.6 


[123] 

Q^qW 

CMS 

8 

19.7 

B2G-12-017 

[124] 

Q^jZ 

CDF 

1.96 

1.055 


[125] 

u ^tz 

CMS 

7 

5 

B2G-12-004 

[121] 


CMS 

7 

1.1 

EXO-11-005 

[126] 


ATLAS 

8 

20.3 

CONF-2015-012 

[127] 

n hn 

CMS 

8 

19.8 

B2G-12-019 

[128] 


CMS 

8 

19.5 

B2G-13-003 

[129] 


CMS 

8 

19.7 

B2G-14-001 

[130] 


CMS 

7 

5 

EXO-11-066 

[131] 

D^bZ 

CMS 

8 

19.8 

B2G-12-019 

[128] 


CMS 

8 

19.5 

B2G-13-003 

[129] 


CMS 

8 

19.6 

B2G-12-021 

[132] 


ATLAS 

8 

20.3 

EXOT-2013-16 

[133] 

D^tW 

CMS 

8 

19.8 

B2G-12-019 

[128] 


CMS 

8 

19.5 

B2G-13-003 

[129] 


CDF 

1.96 

2.7 


[134] 


CMS 

7 

5 

EXO-11-050 

[135] 


CMS 

7 

5 

EXO-11-099 

[136] 

Q^bW 

ATLAS 

7 

4.7 

EXOT-2012-07 

[137] 


ATLAS 

8 

20.3 

CONF-2015-012 

[127] 


CMS 

8 

19.7 

B2G-12-017 

[124] 


ATLAS 

8 

20.3 

EXOT-2013-16 

[133] 

Qhjz tw 

ATLAS 

8 

20.3 

EXOT-2014-17 

[138] 


CMS 

8 

19.6 

B2G-12-012 

[139] 

U ^tH 

CMS 

8 

19.7 

B2G-12-004 

[140] 


Table 2: Experimental analyses included in our numerics for heavy quark partner decay. Q 
stands for any quark partner where the decay in question is allowed by electric charges, 
j stands for a light quark or 6 jet, and q for a light quark jet. 
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• Single production, that is relevant for quarks that have a significant degree of 
compositeness, and in this case dominates at higher masses [26-30,142]; 

• Contributions to the pair production cross section from heavy resonance exchange 
[143,144], 

Taking into account these two effects is beyond the scope of our study, as it cannot be 
implemented efficiently in a fast parameter scan. The bounds we obtain should thus be 
considered conservative. 

Since the experimental analyses typically quote bounds on the pair production cross 
section assuming a 100% branching ratio to the desired final state, we correct for the 
branching ratio BR((5 —)• /) of the quark partner to final state / by multiplying the 
production cross section with 

• BR(Q —)• /)^ in case the experimental analysis requires both partners to decay 

to /; 

• ^1 — (1 — BR((2 —?■ /))^^ in case the experimental analysis requires one or both of 
the partners to decay to /. 

In the M4DCHM5, there are in total 24 heavy charge-2/3 quarks (denoted with U 
in the following), 24 charge-(—1/3) quarks {D), as well as 6 exotic charge-5/3 quarks 
(Q5/3) and 6 charge-(—4/3) quarks (Q_4/3)- The decay modes always involve one SM 
quark or quark resonance plus one kk, Z, Higgs, or vector resonance. In our numerical 
analysis, we can only impose constraints on decays involving SM particles only. This is 
not a strong restriction since the lightest quark partners are always required to decay to 
SM states for kinematic reasons. In table 2, we list all the experimental searches that 
we include in our numerical analysis for the individual decay modes. In this table, Q 
stands for any quark partner where the decay in question is allowed by electric charges, 
j stands for a light quark or b jet, and q for a light quark jet. Note that there are no 
dedicated searches for the Q_ 4 / 3 , but searches of the type Q —)■ {hW,jW,qW) are also 
sensitive to these states. 

An important point concerning the experimental coverage of parameter space is that 
the experiments typically employ a hard pT cut to reduce backgrounds. This maximizes 
the sensitivity to heavy states, but misses out on the low end of the mass spectrum. In 
fact, combining all existing analyses in table 2, we have identified a number of gaps in 
the coverage of quark partner masses. This is illustrated by the plots in figure 2. They 
show the 95% C.L. upper bound on the branching ratio in the decay modes to Vk or Z 
bosons as a function of the quark partner mass. We make the following observations. 

• When kinematically allowed, there is a gap between the LEP bound of 100 GeV 
and the lower end of the Tevatron bounds. This should however not be taken 
seriously as quark partners with mass comparable to the top quark would very 
likely have shown up already. 
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BR(Q ^ qW) BR(Q 5/3 ^ tW) BR(;7 ^ tZ) BR(D ^ tW) 



I I kin. forb. 



300 400 500 600 700 800 900 


mo [GeV] 



mo [GeV] 





mq [GeV] mg [GeV] 


Figure 2: Upper bounds on the branching ratios of quark partners decaying to SM states from 
individual Tevatron and LHC searches. QCD pair production is assumed. Q stands 
for any quark partner where the decay is allowed by the electric charges, q stands for 
a light quark (excluding the b). 
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• More seriously, there are gaps between the upper end of the Tevatron exclusions 


and the lower end of the LHC exclusions. This leaves a window between about 
300 and 350-500 GeV not covered by existing searches. The only exception is 


the mode U —)• tZ, but quark partners around 300 GeV would have a very small 
phase space to decay to tZ, making it plausible that the branching ratio is smaller 
than in the other channels. 

• Bounds on light quark partners are weak, with no existing LHC search for the 
decay mode involving the Z boson. This is problematic since, depending on their 
quantum numbers, some of the light generation partners have very small branching 
ratio into qW as will be demonstrated in section 4.5.1 below. 

Concerning the gaps mentioned in the second item, they could very likely be closed by 
a reanalysis of existing run-1 data (this is also indicated by recasting of some existing 
new physics searches, see e.g. [145]). We call upon the experimental collaborations to 
perform such a reanalysis, given the importance of the partner mass scale for naturalness 
and radiative electroweak symmetry breaking in the models at hand. In our numerical 
analysis, in order not to be biased by these low-mass regions for quark partners, we have 
imposed a hard lower bound of 500 GeV on all quark partner masses, in addition to 
the LHC 7 and 8 TeV searches that are sensitive to higher masses (while the Tevatron 
searches become irrelevant). 

3.3.2. Spin-1 partners 

Spin-1 resonances can be pair-produced in a Drell-Yan like process. If narrow enough, 
they would show up as a peak in the dilepton, dijet, ti, VV, or Vh final state, depending 
on the branching ratios. In the models considered by us, the spin-1 and spin-1/2 res¬ 
onances are strongly coupled to each other. Consequently, if the decay to two fermion 
resonances is kinematically allowed, the resonances become very broad and are not cap¬ 
tured by the experimental analyses anymore^^. Still, we expect that they are sufficiently 
narrow in part of the parameter space, so we include the experimental constraints. 

The hadronic production cross section of a spin-1 resonance can be written as 



(72) 


where qi ^2 = u or d, Cp is the colour multiplicity of the resonance (cp = 1 for electroweak 
resonances and Cp = 8 for pc'), s is the hadronic center-of-mass energy squared and Tipa 
is the parton-parton luminosity of the relevant initial state defined as 



(73) 


^^The presence of a gap between Tevatron and LHC 7 TeV searches has already been noticed in [23]. 
^^They would however still constitute a contribution to the pair production cross section of the fermion 


resonances, see the comment in section 3.3.1 
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Decay 

Experiment 

[TeV] 

Lum. [fb 

Analysis 



ATLAS 

7 

4.7 

EXOT-2012-02 

[146] 

p± ^ w^h 

ATLAS 

8 

20.3 

EXOT-2013-23 

[147] 

CMS 

8 

19.7 

EXO-14-010 

[148] 


ATLAS 

8 

20.3 

EXOT-2013-01 

[149] 

p± ^ 

ATLAS 

8 

20.3 

EXOT-2013-07 

[150] 

ATLAS 

8 

20.3 

EXOT-2013-08 

[151] 


CMS 

8 

19.7 

EXO-12-024 

[152] 

t 

-H 

CMS 

8 

19.5 

B2G-12-010 

[153] 

pO ^ w+W- 

ATLAS 

8 

20.3 

EXOT-2013-01 

[149] 


CMS 

8 

19.7 

EXO-13-009 

[154] 

^ Zh 

ATLAS 

8 

20.3 

EXOT-2013-23 

[147] 

CMS 

8 

19.7 

EXO-13-007 

[155] 

^ £+£- 

ATLAS 

8 

20.3 

EXOT-2012-23 

[156] 

CMS 

8 

20.6 

EXO12061 

[157] 

P^/PG -t ti 

ATLAS 

8 

20.3 

CONE-2015-009 

[158] 

CMS 

8 

19.5 

B2G-12-008 

[159] 


Table 3: Experimental analyses included in our numerics for heavy vector resonance decay. 


For the Tevatron, (72) is valid with the appropriate replacements. 

The M4DCHM5 contains 3 charged and 5 neutral electroweak vector resonances, plus 
the colour-octet gluon resonance. As in the spin-^ case, we can only impose constraints 
on decays to SM states. We include all the analyses listed in table 3. The only relevant 
final state that we have not included is the decay to dijets. The reason is that the 
experimental bounds on the dijet resonance cross section depend on an acceptance factor 
that is model dependent and that is not easy to take into account in a parameter scan. 

In our numerical analysis, we employ a cut of 5% on the relative total width Tp/rup, 
above which all bounds are ignored for an individual resonance. While several of the 
analyses are actually sensitive to broader resonances, it is not possible to include this 
dependence in a parameter scan. As in the case of fermion resonances, our bounds 
should thus be considered as conservative. 

3.3.3. LHC excesses 

Interestingly, several of the searches for spin-1 resonances we include as constraints 
contain an excess of events around 2 TeV. The most significant deviation is in the 
ATLAS search for —)• WZ, corresponding to a local significance of 3.4cj, but an 

excess around the same mass appears also in the corresponding CMS search, and, to a 
lesser extent, in searches for WW oi Wh final states. This is particularly interesting in 
the context of composite Higgs models as the resonances associated to SU (2 )i (denoted 
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PLfj, in the gauge eigenstate basis in section 2) form a triplet of a charged and a neutral 
vector that are almost degenerate, have approximately equal branching ratios to WZ, 
WW, Wh, and Zh final states, and can have a production cross section in the right 
ballpark to explain these excesses [36,37,160-164], 

4. Numerical analysis and predictions 

This section contains the main results of our paper. After describing the numerical 
analysis procedure in section 4.1, we will discuss fine-tuning in all models in section 4.2.2, 
the numerical results of signals in indirect searches in U{2)f^Q in section 4.3, for indirect 
searches in U{2)^q and in section 4.4, and for direct LHC searches in all models 

in section 4.5. 

4.1. Strategy 

4.1.1. Scanning procedure 

Our aim is to sample the parameter space of the M4DCHM with four different flavour 
structures while satisfying all the experimental constraints discussed in section 3. This 
is particularly challenging because partial compositeness implies that all SM masses and 
couplings are relatively complicated functions of the model parameters and the additional 
requirement of correct radiative EWSB is even harder to control analytically. To cope 
with these challenges, we have improved a method first employed in [46]. We construct 
a function depending on all the theoretical predictions for all constraints discussed 
in section 3 as well as the corresponding experimental measurements. We then proceed 
in four steps. 

1. We randomly generate a set of input parameters that only fulfills the most rudi¬ 
mentary consistency conditions (e.g. composite gauge couplings being greater than 
1, effective potential possesses a minimum away from zero). 

2. We use brute-force numerical minimization (with NLopt [165]) to “burn-in” into a 
region of parameter space not too far from viability. 

3. We use a Markov Chain Monte Carlo (MCMC, with pypmc [166]) to sample the 
viable parameter space. 

4. We filter the Markov chains so that only points remain where each individual 
constraint is satisfied at the 3a level. 

After the burn-in phase, the generation of viable parameter points is very efficient, as 
the MCMC is adaptive and has an acceptance rate around 23%. The downside of the 
method is that adjacent points have high autocorrelation, implying that very long chains 
are needed to obtain a reasonable coverage of the parameter space. Furthermore, the 
parameter space can contain several disconnected minima. For these reasons, we run a 
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large number of chains (around 500 for each model) starting at different (random) initial 
points. 

We stress that we do not interpret the outcome of the Markov Chain statistically, 
in the sense of a posterior probability distribution for the model parameters. Apart 
from the problem of reaching sufficient coverage of the parameter space, this would be 
problematic due to dependence on the choice of priors. Instead, we use the MCMC 
algorithm as a shortcut to generate a sufficient number of valid points. In our final 
sample, these points approximately follow a normal distribution peaked around 30-40, 
for 48 individual contributions to the We also find that model-independently for 
nearly all points only < 5 individual constraints are violated by more than 2a at a time. 
Consequently, the hard 2>a cut only removes a small fraction of extreme points. The 
largest deviations are typically found for m*, BR(i? —)• Xg'y), the inclusive values of Vub 
and Vcb, and for the Higgs signal strengths. 

4.1.2. Model parameters 

Below, we specify the model parameters and any relations we have imposed among them 
in our scan. 

• /) fi, fx, fc- 

We have imposed 1 < /i// < Vs, where the lower bound corresponds to a de¬ 
coupling of the axial resonances and the upper bound is motivated by the partial 
unitarization of Goldstone boson scattering [13]. We have not assumed /i, fx, 
and fc to be degenerate, but we have restricted them to be within a factor of two, 
i-e. 2 ^ /x,g//i < 2, to prevent the fit from completely decoupling one of the 
resonances. 


• 9p, 9x, 9G- 

We have varied these couplings completely independently, only imposing 1 < 
9p,x,G < 47r to have a strong but semi-perturbative coupling (in the case of the 
gluon resonance, we also imposed gc > 9so)- We further imposed /i,x,G 9p,x,G /< 
47r/, to not have resonance masses above the (naive) cutoff of the theory. 

• 9o, 90, 9s0- 

These parameters are fixed by the known gauge couplings once the composite gauge 
couplings are specified. 

• TTiQ, rUg, oiyq, txlyq + Yq where Q = [/ or Zl, in the case of U{2)^ different for the 
first two and the third generation. 

In our numerical analysis described above, we have treated the logarithms of these 
parameters as scan parameters, in order not to get a bias towards heavier masses. 
The only relation (apart from the ones forced by the flavour symmetries) we have 
imposed is that all these parameters are < dvr/. Note that this implies that the 
above parameters can only take positive values in our scan, but this can always be 
arranged by a suitable choice of phases for the fermion fields. 
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• Quantities parametrizing the composite-elementary mixings, see appendix C for 
details. 

Again, we have scanned dimensionful parameters logarithmically and require them 
to be < 47r/, but otherwise we impose no restrictions (many relations among these 
parameters are fixed by the requirement to have the correct quark masses and 
CKM mixing). 

With these assumptions, the total number of real parameters or phases is 44 for U{2)'^q 
and U{2)\q and 30 for and To compare these parameters to the SM, 

it should be noted that we do not treat lepton masses as free parameters, have massless 
neutrinos and set the QCD 9 term to zero, but the Higgs mass and VEV are predictions 
rather than inputs. 


4.2. General fit results and fine-tuning 
4.2.1. Failure of f7(3)Lc 

In the case of U{3)f^Q, our scans have not been able to find a single viable parameter 
point, even for a large number of chains. This is not surprising as already a qualitative 
analysis of the relevant constraints on [/(3 )lc [18] has revealed that there are extremely 
strong constraints on the model from electroweak precision test and CKM unitarity. It 
seems to be impossible to reconcile these constraints with the need for correct radiative 
EWSB. We will thus not consider U{3)j^q any further. 


4.2.2. Fine-tuning 


Before discussing predictions for physical observables, let us address the degree to which 
the viable model points we have found can address the electroweak hierarchy problem. 
To this end, we have computed the Barbieri-Giudice fine-tuning measure [167] 


Abg = max 

i 


din mz 
d In Xi 


(74) 


that quantifies the sensitivity of the weak scale to variations in the model parameters 
Xi- In composite Higgs models, Abg can be obtained directly from derivatives of the 
potential [168] . Still, given the large parameter space, the computation turns out to be 
more time-consuming than for the physical observables, so we have computed Abg only 
for a subset (2%) of all our points. The results for the three viable models are shown 
in fig. 3.^^^ Not surprisingly, the lowest Abg is obtained for low / and a sub-percent 
fine-tuning is possible in all models as long as / < 1 TeV. This is compatible with 
earlier analyses in similar models [168,169]. The points with the lowest tuning measure, 
highlighted by stars in the plot, have 

. A^a^ = 33 for U{2)lc, 

^^Note that the individual “speckles” - visible in many of the scatter plots - correspond to different 

Markov chains. 
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. = 55 for [/(2)|c, 

. = 73 for t/(3)|c. 

Two comments are in order here. First, we reiterate that we define viability for a point 
as fulhlling all individual constraints at 3a. Since we assume a 5% relative uncertainty 
on ruf and mh (see section 3.1), the known tendency of the model to have a too heavy 
Higgs and a too light top means that the points with lowest tuning typically have the 
Higgs and top mass 15% away from their central values. Second, we stress that there are 
variants of the model considered by us that have lower hne-tuning since the M4DCHM5 
suffers from a “double tuning” by an accidental enhancement of the Higgs mass due to 
the structure of the potential, see [168] for a discussion and alternatives. 

To get a better understanding of the tuning in the Higgs potential, we adopted an 
expansion of the potential in terms of Sh (as also used e.g. in [13]), 

Ves--lf^sl + /3f.sl (75) 

where we have explicitly defined the 7 and /3 parameters as dimensionless by factoring 
out their typical scale and we have neglected terms of 0(s®). In this notation, the 
Goldstone VEV and the Higgs mass are given as 

=\/i' - i) f"- 

The requirement to fulhll EWPT’s (and therefore to have a not too large s)j) together 
with the hierarchy rrih / forces 7 to take a rather small value. As already mentioned 
in [13,168], this requires a cancellation between the fermion and gauge contributions to 
a large degree. 

In our framework, we can extract the 7 and /3 parameters for each contributing field 
by simply htting the numerical values of (23) to the parametrization (75). We indeed 
find that the fit prefers highly correlated gauge and fermion contributions that are large 
individually but almost completely compensate each other. We also find large cancella¬ 
tions between the up- and down-quark sector and also between individual contributions 
in each sector. 

4.3. Left-handed compositeness: indirect searches 

As discussed in section 4.2.1, we have not found any viable points for U{3)f^Q. We will 
thus restrict ourselves to U{2)\^q in this section. 

4.3.1. Light quark compositeness 

Compositeness of the first two generation quarks is mainly constrained by 

• Eirst-row CKM unitarity, see section 3.1.2; 

• The hadronic Z width, see section 3.2.2; 
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Figure 3: Barbieri-Giudice fine-tuning measure vs. / for the three viable models for a thinned- 
out sample of all our viable parameter points. The stars show the points with lowest 
fine-tuning measure for each model. 



Figure 4: Observables sensitive to light quark compositeness in 17(2)^^. Left: first-row effective 
CKM elements. The black line corresponds to the SM limit of a unitary CKM matrix. 
Right: hadronic Z width (normalized to Z ^ ee) vs. the pp —>■ jj cross section in 
the rapidity bin described in appendix D, normalized to the 95% C.L. limit extracted 
from the ATLAS analysis. The black line corresponds to the central value of the SM 
prediction. 
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• The dijet angular distribution at LHC, see section 3.2.6. 

The predictions for these quantities are shown in figure 4. The left-hand plot shows 
the effective CKM elements |14s| vs. \Vud\ and demonstrates that large deviations from 
the SM relation |14sp -|- \ ^ 1 (shown as a black line) are possible. The solid gray 

lines show the current experimental 2 cj bounds. At 3a, the points stop because of our 
procedure described above. On the right, we show the predictions for the hadronic Z 
width as defined in (37) (the SM central value is shown as a black line) vs. the pp —)• jj 
cross section in the bin described in appendix D, normalized to the 95% C.L. limit 
extracted from the ATLAS analysis. This plot shows that sizable effects are possible 
in these observables as well, but almost all points lie within the 2a region for both 
observables (shown as solid gray lines), demonstrating that CKM unitarity is by far the 
strongest constraint on light quark compositeness in ?7(2 )lq at present. 

In these plots (as in almost all the plots of this section), on top of all the viable 
points in blue, we show points that have a fine-tuning measure Abg < 100 in yellow. 
The rationale is to demonstrate in which part of the viable space for the observables 
in question the most “natural” parameter points lie. We warn the reader however that 
these points do not correspond to all viable points with Abg < 100 - as mentioned in 
section 4.2.2, we have only computed Abg for a subset of the points. One should also 
keep in mind that, simply due to their smaller number, these points typically cluster in 
the region with the highest point density. 

4.3.2. Higgs production and decay 

The left-hand plot in figure 5 shows the signal strengths of the Higgs produced in gluon 
fusion and decaying to ZZ (which equals the one to WW due to custodial symmetry), 
77 , and bb. Most of the points lie on the curves that are expected from analytical 
considerations of coupling modifications (see e.g. [170]). This is even true for h —)• 77 
and gg —)• h, since in pNGB Higgs models, the loop contribution of heavy resonances to 
these processes almost entirely cancels with the coupling modification of the top quark, 
leaving the Higgs non-linearities as the dominant effect [22]. 

However the plot also shows deviations from these relations. These can be understood 
to be caused by light quark compositeness, spoiling the above mentioned cancellation 
[171]. In this way, the signal strength can be closer to (or further away from) their SM 
values than naively expected for small /. Nevertheless, we find this effect to be mild, 
given the strong constraints on light quark compositeness discussed in the previous 
section. 

Concerning the h ^ bb signal strength, we note that the figure shows the signal 
strength in the case of gluon fusion production, while the experimental bounds are 
currently based on the associated production with vector bosons, that we do not include 
in our analysis. 


^®The fact that almost all points lie on the same curve and only a few individual Markov chains have 
strayed away from it explains the frayed appearance of the plot. 
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Figure 5: Left: Higgs signal strength for gluon fusion production and decay to final states ZZ 
(equal to WW by custodial symmetry), 77 , and bb in U{2)f^Q. The SM corresponds 
to /i = 1, shown as a horizontal line. Right: Oblique parameters S and T in {7(2 )lc, 
defined to be 0 in the SM. 


4.3.3. Oblique parameters 

The right-hand plot in figure 5 shows the predictions for the S and T parameters. 
We show the region allowed by experiment at 2cr as a gray dashed ellipse, while the 
gray solid ellipse takes into account also the additional theory uncertainty discussed in 
section 3.2.1. The tree-level contribution to S is strictly positive, while the fermionic 
loop contribution to T can have either sign, but is preferred to be positive by experiment 
for positive values of S and indeed large positive contributions are possible for our choice 
of fermion representations, which is important as it helps to alleviate the bound from S. 

4.3.4. Meson-antimeson mixing 

Figure 6 shows the predictions for AT = 2 observables in U{2)y^q. In the left-hand col¬ 
umn, we directly show the correlation between observables. In this case, it is important 
to notice that the CKM parameters themselves are varied in our fit and are not fixed to 
their SM central values (as is often done in parameter scans of, e.g., SUSY models). As 
a consequence, any correlations among (tree-level) NP contributions are washed out by 
the spread in the allowed values for the CKM parameters. The dashed gray lines in these 
plots show the allowed regions with merely the experimental 2a uncertainties, while the 
solid gray contours take into account additionally the (correlated) theory uncertainties 
at 2a. As discussed in section 3.2.4, we only need to account for the non-CKM (i.e., 
lattice) theory uncertainties, as for a given point, the CKM parameters are predictions. 
While these plots are more directly related to the experimental measurements, the 
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Figure 6: AF = 2 observables in [/(2 )lc 
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variation of CKM parameters obscures the relation to the NP contributions. This is 
why in the right-hand column of figure 6 we show the ratios (or phase differences) of 
the observables and the SM VP-loop contribution for each parameter point. In this way, 
the correlations valid at leading order in the C/(2)^ spurion expansion, AM^/AMJ^ = 
AMs/AMf^ and 4>s — 4>f^ = 4>d — (where S^Ks = siiKpd), shown by solid black 
lines in the plots, become apparent. In the bottom right plot, the black line corresponds 
to the MFV limit of equal relative modification in the and mixing amplitudes, 
while no such correlation is expected in U{2)^. We make the following observations. 

• The mass differences in the Bd and Bg systems can receive corrections up to -|-60%, 
but negative NP contributions are disfavoured. This can be understood from the 
fact that the tree-level Wilson coefficient of Qy^i, cf. (44), involves the square of 
a coupling that carries a small phase. 

• The Bg mixing phase can saturate the experimental lower bound, but positive 
values for siiKpg are only predicted for a small number of points. This is due to 
the correlation with c/d and the preference for a negative NP contribution to the 
latter, that is also visible in global CKM hts [80,172,173]. 

• Both for the mass differences and for the phases in the Bd and Bg systems, the 
leading-order U{2)^ correlations are broken for a significant fraction of the points, 
seen as a deviation from the black diagonal lines. This is due to non-negligible 
contributions from left-right operators. We have identified two reasons for why 
these effects are larger than expected from a general EFT analysis [61]. 

1. The Wilson coefficients of these operators are RG-enhanced; 

2. Due to partial compositeness and the possibility to have a hierarchy even 
among the (diagonal) left-handed composite-elementary mixings, the spurion 
hierarchies in the right-handed mixings can be milder than the Yukawa hier¬ 
archies, effectively enhancing subleading terms in the spurion expansion. 

A similar effect has already been noted in the context of the MSSM with a U{2)^ 
symmetry [172] (where it was mostly due to an accidental enhancement of a loop 
function) and we find the effect to be even more pronounced in the composite 
Higgs case. We stress nevertheless that the majority of parameter points does 
fulhll the U{2)^ relations to a good precision, corresponding to a large density of 
points around the black lines in the plots. 

• The relative modification of ex compared to the SM is always equal'^® to or smaller 
than the relative modification of AMd- This confirms the general expectation for 
Ui2)lc in [18]. 

So far, we have not discussed D^-D^ mixing. On the one hand, the system is 
plagued by large theoretical uncertainties due to poorly known long-distance contribu¬ 
tions; on the other, the effects in U{2)^ models are expected to be small on general 

^®With a small correction factor stemming from the SM charm contribution. 
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grounds [61]. To investigate whether this expectation is correct, we have computed the 
tree-level NP contribution to the mixing amplitude in U{2)f^Q. Since the SM contri¬ 
bution is expected to be real to a good accuracy, the most promising NP effect would be 
a CP violating one. A global fit to data from the D system [174] allows to directly con¬ 
strain the absolute value and the phase of the mixing amplitude. At 2a, this constrains 
the imaginary part of the mixing amplitude to be 

— 0.5 ns“^ > ImM;j^ < 1.6 ns“^. (77) 

Numerically, we have found that the NP contributions to ImM^^ are always negative 
in C/(2 )lc and can reach at most —0.5 ns“^. We conclude that CP violation in 
mixing is currently not a relevant constraint on the model, but future improvements of 
the bound by factors of a few would start to cut into its parameter space. We have 
further found that the NP contributions to ImMj^ are strongly anticorrelated with 
eK'- sizable NP contributions to the former never occur simultaneously with sizable NP 
contributions to the latter. However, the NP contributions to both observables can be 
small simultaneously. 

4.3.5. Rare B decays 

The Wilson coefficient Cj^ of the electromagnetic dipole operator, cf. sec. 3.2.5, receives 
NP contributions, but only to the extent that is allowed by the strong constraint from 
the branching ratio of H —)• Xs'j- We find these contributions to be aligned in phase with 
the SM to a high degree, such that CP violating effects, e.g. in the direct CP asymmetry 
in H —)• K*^, are expected to be small. Contributions to the chirality-flipped coefficient 
are small by U{2)^ symmetry. 

The most interesting effects in rare B decays stem from the tree-level contributions to 
the semi-leptonic Wilson coefficients Cg^ and Cig. As discussed in sec. 3.2.5, there are 
Z-mediated and resonance-mediated effects that dominantly contribute to Cj'q, but also 
resonance-mediated effects that contribute only to Cg^. We remind the reader that in our 
numerical analysis, the only observable sensitive to these Wilson coefficients that we have 
imposed as a constraint is the branching ratio of Bg —)• //"’"/U”, essentially limiting the 
absolute value of C^g. All points passing this constraint are shown in the left-hand plot 
of fig. 7. We observe that large NP effects in C^g - saturating the experimental bound 
on Bg —)■ - are possible, but also sizable effects in Cg^. Interestingly, the largest 

effects allowed in Cg^ correspond to a negative sign that is preferred by the anomalies 
in B —)■ angular observables discussed in sec. 3.2.5. The gray ellipse in figure 7 

left corresponds to the 2a preferred region in a global fit to 6 —)• sfi'^ observables [115], 
which shows a clear tension with the SM point (0,0). The figure clearly shows that if 
these tensions are due to NP, the M4DCHM5-t/(2)Lc can explain them. This is also 
demonstrated by the right-hand plot in fig. 7, which shows the predictions for two of 
the observables that currently show the biggest tensions with the SM, namely the low-g^ 
branching ratio of Bg —)• (/>/i"*■//“ and the angular observable P 5 in R —)• K*. In this 
plot, the black star shows the central value of the SM predictions (taken from [105,115]), 
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Figure 7: Left: new physics contributions to the Wilson coefficients Cg^ and C\q in 

Right: predictions for the angular observable P 5 in R —>■ and the branching 

ratio of Bg —>■ (pfJL^, both in the low-g^ bin from 1 to 6 GeV^. The star corresponds 
to the central values of the SM predictions. 


the gray dashed line the values allowed at 2(t by experiment, and the gray solid lines the 
2a allowed values taking into account also the theoretical uncertainties. 

We have found that all of the points that have Cg® < —0.5 - and could thus account 
for the tensions in angular observables and branching ratios - correspond to a small 
value (between 1 and 2) of the composite coupling gx and a correspondingly small mass 
(below 1 TeV) of the mass eigenstate that is dominantly the px resonance. This can be 
understood from the discussion in section 3.2.5: in the limit gx ^ Qp, the “KK photon”- 
like state is dominantly the field px- In addition, since gx is not much larger than 
the elementary gauge coupling, the parametric suppression of the resonance-mediated 
contribution is lifted. It is also important to notice that this linear combination of gauge 
fields does not contribute to the S parameter and thus is allowed to be lighter than the 
other vector resonances. Sizable contributions to Bg-Bg mixing are also generated by 
the exchange of the light resonance, but we find that the shift in AMg is below 20% 
relative to the SM. We have also computed the LHC production cross section and decay 
branching ratios of the light resonance for the points with sizable NP effects in Cg. For 
most of the points, the dominant decay mode is tt and the cross-section is just below 
the ATLAS and CMS searches for resonances in this mode that we have imposed as a 
constraint in our scan (see section 3.3.2). The width of the resonance is small enough to 
show up in a “bump hunt”. Prospects for vector resonances will be discussed in more 
detail in section 4.5.3 below. 

In the left-hand plot of fig. 8, we finally show the predictions for the correlation 
between the branching ratios of Bg —)• p~ and B^ —)• which is fixed by U{2)^ 
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Figure 8: Left: predictions for the branching ratios of Bg —)■ and —>■ in C/(2)lq. 

Right: predictions for the deviation of the effective CKM element Vtb from 1 vs. the 
branching ratio of the FCNC top decay t —>■ cZ in U{2)f^Q. 


to be equal relative to the respective SM predictions (but is again slightly washed out 
by the variation in CKM elements). The current 3cr upper bound on Bg —?■ can 

be saturated, but also a significant suppression can occur. This is in contrast to, e.g., 
the Littlest Higgs model with T-parity, where this branching ratio can only be enhanced 
with respect to the SM [175]. 

4.3.6. Top decays 

A significant degree of compositeness of the left-handed top quark can lead to a reduction 
of the single top production cross section at LHC, corresponding to a reduced value for 
the effective CKM element Vth as discussed in section 3.1.2. In addition, in this case 
there can be sizable flavour-changing couplings of the top quark to the Z boson, since 
the left-handed couplings are not custodially protected, in contrast to the right-handed 
ones. 

These two effects manifest themselves in a correlation between the deviation of Vth 
from 1 and the branching ratio of the FCNC top decay t —)• cZ as shown in the right- 
hand plot of figure 8. Both effects are quite moderate after imposing all the bounds. The 
deviation in Vth is always within the current 2cr experimental constraint and percent-level 
experimental accuracy will be necessary to find a significant deviation. The branching 
ratio of the FCNC top decay can reach at most 10“^, which will be challenging to see 
at the LHC [20,176]. 
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Figure 9: Hadronic Z width (normalized to Z —>■ ee) vs. the pp —>■ jj cross section in the rapidity 
bin described in the text, normalized to the 95% C.L. limit extracted from the ATLAS 
analysis. The black line corresponds to the central value of the SM prediction. Left: 
^(2 )rc> right: 

4.3.7. Other processes 

So far, we have not discussed rare K decays. While these processes are important 
constraints on many NP models, we find the effects in U{2)f^Q to be rather small. For 
instance, the branching ratios of K~^ —)• and Kp —)• are modified by at most 

±20% with respect to the SM (and are perfectly correlated due to U{2)^). The short- 
distance contribution to the branching ratio of Kp —?■ pL~ is always below 2 x 10“®. 

4.4. Right-handed compositeness: indirect searches 

In contrast to C/(3 )lc, we do find a viable parameter space for the U{2>)\q model. Since 
C/(3)^c is a limiting case of the more general U{2)^q (the limit in which the composite 
sector mass parameters and the right-handed composite-elementary mixings for the first 
two and the third generation coincide), it is natural to discuss them together. We will 
proceed as in the case of left-handed compositeness in section 4.3. 

4.4.1. Light quark compositeness 

In the case of right-handed compositeness, it is typically the right-handed light quarks 
that can carry a sizable degree of compositeness. Consequently, in contrast to left-handed 
compositeness, first-row CKM unitarity does not constitute a relevant constraint and the 
main constraint is given by the hadronic Z width and the dijet angular distribution. The 
predictions for these quantities are shown in figure 9 that is the analogue of figure 4 right. 
We make the following observations. 
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Figure 10: Higgs signal strength for gluon fusion production and decay to final states ZZ (equal 
to WW by custodial symmetry), 77 , and bb in U{2)^(. (left) and U{3)^q (right). 
The SM corresponds to /r = 1, shown as a horizontal line. 


• In both models, Rh is within the 2cj bounds for almost all the points. 

• In both models^^, large effects relative to the experimental constraints are ob¬ 
tained in the dijet angular distribution. This is the strongest bound on light-quark 
compositeness in the right-handed compositeness models. 

• In the case of U{3)^q, we see that there is even a /ower bound on the modification 
of the dijet angular distribution. Improved experimental measurements in the 
future could help to disfavour this scenario. 

4.4.2. Higgs production and decay 

Figure 10 shows the Higgs signal strengths for right-handed compositeness in analogy 
to figure 5. As discussed in section 4.3.2, the leading dependence on / is modified by 
light quark compositeness, which is more pronounced in U{3)^q due to the requirement 
to have a large degree of compositeness for all right-handed up-type quarks. 
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Figure 11: AF = 2 observables in U(2)^q 


42 
























































4.4.3. Meson-antimeson mixing 

Figure 11 shows the predictions for AF = 2 observables in U (2)^q, in analogy to figure 6. 
We first point out the similar features, 

• Sizable enhancements at the level of 60% with respect to the SM are possible in 
the Bd and Bg mass differences, but a suppression is strongly disfavoured. 

• The leading-order U{2)^ relation between AM^ and AMg (shown as a black line) 
is violated by LR operators. 

But there are also important differences between 1/(2)^^ and U{2)^q. 

• There is no new phase in B^ mixing, as was already pointed out in [18]. 

• In Bg mixing, on the other hand, there can be a new phase roughly at the level 
of the current experimental uncertainties. This phase stems from the subleading 
terms in the spurion expansion and thus violates the leading order U{2)^ relation 
(implying equal phase shifts in B^ and Bg mixing). 

• The enhancement of ex relative to the SM is always larger than the one in AM^- 
This is the opposite of what happened in ?7(2)lc, where the relative enhancement 
was always smaller in ex, cf. figure 6 bottom-right. In the future, this could serve 
as a way to distinguish the two models based on AF = 2 observables alone. 

In U{3)^q, we only directly show the observables normalized to their SM values in 
figure 12. In this case, the MFV relations, shown by black lines, are fulfilled exactly and 
there is no new phase, neither in B^ nor in Bg mixing. 

Concerning D^-D^ mixing, in [/(2)^q, similarly to U{2)f^Q discussed at the end of 
section 4.3.4, the NP effects are quite small and we find that the imaginary part of the 
mixing amplitude is always between —0.4 and -1-0.2 ns“^, which is not relevant at the 
current experimental precision, but will become relevant when the experimental bound 
improves by an order of magnitude. In U{3)^q, there is no new phase and thus no NP 
contribution to the imaginary part of the mixing amplitude. 

4.4.4. Rare B decays 

In U{2)^fj, similarly to the case of left-handed compositeness discussed in section 4.3.5, 
the largest contribution to the Wilson coefficients of the semi-leptonic b —)• sii transition 
occurs in the Wilson coefficient C^q, but there are also contributions to the Wilson 
coefficient as shown in figure 13 left. In this case, we only find a small number 

^^The alert reader may have noticed that in P(3)pic, many points saturate the experimental upper 
bound, while in U{2)^q, this does not seem to be the case, even though we have stated that U{3)^q 
is a subset of the U{2)^q model. The reason is a volume effect in the high-dimensional parameter 
space; one would need a huge number of points in the U{2)^q model to get a reasonable coverage of 
the subspace corresponding to U{3)^q. This effect is visible in many of the plots in this section and 
justifies the separate analysis of U{3)^c- 
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Figure 13: Left: new physics contributions to the Wilson coefficients and CJo in 17(2 )|^q. 

Right: predictions for the branching ratios of Bg —>■ in 
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Figure 14: Predictions for the angular observable P 5 in i? —)■ K* and the branching ratio 

of Bg —1 both in the low-g^ bin from 1 to 6 GeV^ in t/(2)|^Q (left) and 

G(3)rc (I'ight). The star corresponds to the central values of the SM predictions. 


of points with sizably negative Cg^ that populate the region preferred by a global fit 
to 6 —>■ sU data, indicated by a gray ellipse. These points then predict a significant 
snppression in absolute value of the angular observable P 5 in B —)• at low q'^, 

see the points around Pg ~ —0.1 in figure 14 left. A distinguishing feature compared to 
U ( 2 )lc is that the contributions to C'J'q are almost always positive, implying that the 
branching ratio of Bg —)• is almost always suppressed, as shown in fig. 13 right. 

In U{3)^q, the contributions to CJ’q are forbidden by an interplay between custo¬ 
dial protection and the flavour structure as discussed in section 3.2.5. However, the 
resonance-mediated contributions to Cg^ are still present and we find viable points in 
the range —1.7 < Cg^ < 0.9. Consequently, also U{3)^q can explain the anomalies 
in 5 — 7 - sii angular observables and branching ratios. For the observable P 5 and the 
branching ratio of Bg —t- , this is illustrated in figure 14 right. Finally, we remind 

the reader again the U{3)^q is actually a limiting subset of P( 2 )^q, so the fact that in 
P(3)^c there are much more points with sizable NP effects in Cg^ compared to U(2)^(j 
is simply a statistical effect since the 17(3)^^ parameter space is more restricted. 

As in l7(2)f^Q discussed in section 4.3.5, the solution of the B physics anomalies by a 
negative NP contribution to Cg^ implies the presence of a light, narrow neutral vector 
resonance below about 1 TeV. In U{2)^q, the dominant decay mode of this resonance is 
tt or two light quark jets, while in U{3)^q the dominant decay mode is always dijets. 

4.4.5. Other processes 

We have not discussed the oblique parameters as the predictions in both models with 
right-handed compositeness are analogous to the effects in U{2)j^q shown in fig. 5 right, 
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so the same comments as in section 4.3.3 apply. 

In rare K decays, the effects both in ?7(2 )^q and in U{3)^q are even smaller than in 
C/(2 )lc discussed in section 4.3.7. 

In contrast to left-handed compositeness, the branching ratio of the FCNC top decay 
t —)■ cZ is always below 10 “® in 1 /( 2 )^^ and even below 10 “^ in U{3)^q and thus 
negligible. 

4.5. Direct searches in left- and right-handed compositeness 
4.5.1. Prospects for quark partner searches 

The direct bounds on quark partner masses discussed in section 3.3.1 are among the 
most important constraints in our analysis. It is thus clear that future searches for 
quark partners will be instrumental in probing these models. Since in our numerical 
analysis, the lightest vector resonances are always found to be heavier than the lightest 
quark partners, which is due to electroweak precision tests and the other indirect bounds 
discussed in section 3.2, the lightest quark partners always decay to SM states. To judge 
which of the search channels will be most promising at run 2 of the LHC, let us first 
discuss the dominant decay modes of the lightest resonances. 

Exotic charge quarks The charge-5/3 and charge-(—4/3) quarks always decay to 
a W boson and a SM quark or quark resonance. In 1/(2)^^, we find that there is a 
significant number of points where the Q 5/3 can decay to both Wt and Wq {q = u,c) 
with significant branching ratio, and similarly for the Q-^/^ decaying into final states 
with bottoms vs. light quarks. In the right-handed compositeness models, we find in 
contrast that for any given exotic quark partner, only the decay to 3rd generation or 
the one to light quarks is relevant. This can be understood from the fact that the decay 
of the exotic charge quarks always involves right-handed composite-elementary mixings, 
and these are flavour-diagonal in right-handed compositeness, but involve flavour mixing 
in left-handed compositeness. In figure 15, we show the predictions for the branching 
ratios as a function of the mass for the exotic charged quark partners for a subset of 
all viable points in all three models. An interesting feature of these plots is that there 
is a significant number of points with branching ratios different from zero or one in a 
given channel. Apart from flavour mixing, this is due to the competition with decays 
involving a fermion resonance in the final state. For heavier masses, the branching ratios 
into SM-only states decrease, as can be seen from the plots as well. 

Up- and down-type quark partners When decaying to SM states, these quark 
partners can decay to a W, Z, or h plus a SM quark. In figures 16 and 17, we show 
predictions for the masses and branching ratios in the most important channels for quark 
partners in the three viable models. 

In summary, the plots show that searches for pair-produced quark partners, both with 
exotic and with SM-like charges, are very promising, with masses and branching ratios 
just above what LHC has excluded in run 1 being viable in all models. 
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Figure 15: Branching ratios of charge-(—4/3) (first row) and charge-5/3 (second row) quarks 
to light (left) and third-generation (right) quarks as function of their mass for all 
three models. The coloured regions are the same as in fig. 2. 


4.5.2. LHC excesses 

As discussed in section 3.3.3, several excesses with significances up to around 3 cj have 
been observed by ATLAS and CMS in resonance searches in Wh, WZ, and WW final 
states around a resonance mass of 2 TeV. To investigate whether the models studied by 
us could account for these anomalies, we have computed the production cross sections 
of charged and neutral electroweak vector resonances times the branching ratios to the 
relevant final states. In figures 18 and 19, we show these predictions in the relevant 
mass region for all three viable models, compared to the expected (dashed) and observed 
(solid) limits in some of the relevant ATLAS and CMS searches (for a total list of searches 
included, see section 3.3.2). In these plots, to be conservative we only show points where 
the decaying resonance has a narrow width, namely T/m < 0.05, because, as discussed 
in section 3.3.2, we have not imposed any LHC constraints on broader resonances. We 
note however that there are a significant number of more points in the same region 
where the width is slightly larger than 5%. But even with this strong condition, we do 
find points in all three models where there are resonances with mass around 2 TeV and 
with cross sections of the order of 5 fb in the case of —>• W^h and —)■ W^Z^, 
which is the right ballpark to explain the excesses (see e.g. [160-164]). In the case of 
pO —)• , the predicted cross section is roughly a factor of two smaller due to the 

PDF suppression, but this agrees at least qualitatively with the less pronounced excess 
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Figure 16: Branching ratios of up-type quark partners to final states involving light (left) and 
third generation (right) quarks as function of their mass for all three models. The 
coloured regions are the same as in fig. 2. 
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Figure 17: Branching ratios of down-type quark partners to final states involving light (left) 
and third generation (right) quarks as function of their mass for all three models. 
The coloured regions are the same as in fig. 2. 
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in the CMS analysis, as seen in the upper plot of figure 19. 

Interestingly, a slight excess around 2 TeV has also been observed in a CMS dilepton 
resonance search [157]. Our predictions for this channel are shown in the lower plot of 
figure 19. Also here, we find a significant number of points with cross section times 
branching ratio of the order of 0.1 fb, which could account for this excess. Also in this 
plot, we are only showing resonances with a narrow width. This is also why there are 
few points in the region of interest for the U{3)^q model. In this model, the electroweak 
resonances are typically broader than 5% due to the stronger coupling to light quarks 
compared to the U{2)^ models.^® 

Given figures 18 and 19, the question arises whether the points explaining the excesses 
in the individual plots are actually the same points, i.e. the question whether the models 
can explain all excesses simultaneously. For the hnal states involving bosons, this is 
obviously the case as the branching ratios are sizable only for the composite SU{2)l 
triplet, for which the branching ratios into WW, WZ, and Wh final states are expected 
to be the same (see section 3.3.3). For the diboson vs. dilepton final states, this is not 
obvious, so in figure 20, we compare the cross sections times branching ratios of neutral 
vector resonances decaying to dileptons vs. WW in all cases where the mass is between 
1.7 and 2.2 TeV and the F/m is at most 5%. We observe that the points with production 
cross section times branching ratio into WW of order 1 fb typically lead to a signal in 
dileptons that is one to three orders of magnitude smaller. Comparing this to figure 19, 
we conclude that if the excesses in diboson final states are due to composite resonances, 
the excess in dileptons could be explained as well, but could also be absent. 

4.5.3. Prospects for vector resonance searches 

The discussion in the previous section has already shown that the diboson and dilepton 
final states are promising channels to look for vector resonances in the models studied 
by us. It should however be stressed that the vector resonances are not required to be 
light enough to be probed at LHC, even at y/s = 13 TeV. In all three models, we have 
found viable points with moderate hne-tuning where all vector resonances are heavier 
than 6 TeV. In the following, we discuss the most promising search channels for the 
vector resonances if they are light enough. 

Gluon resonance pc can only decay to fermion pairs and usually has the largest 
branching fraction into quark partners because it couples to them through the strong 
coupling go. In that case, the most promising experimental strategy is to look for the 
quark partner pair decaying to SM particles [143,144]. If the decay to quark partners 
is kinematically disfavoured or forbidden, pc can also decay to SM quark pairs. In the 
U (2)^ models, we find that the decays to ti or bb can be up to 50% and to light quarks 
up to 30% (summing over the four light quark flavours). In 17(3)^^, due to the large 
degree of compositeness of light right-handed quarks, the dominant SM decay mode are 


'^^This does not mean that this model cannot explain the excesses, but a detailed analysis of the impact 
of broad resonances is beyond the scope of our study. 
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Figure 18: Predictions for the production cross sections times branching ratios of charged elec- 
troweak vector resonances decaying to Wh or ZW final states in all three models. 
Only points with narrow resonances (P/to < 0.05) are shown. The dashed and solid 
curves show the expected and observed 95% C.L. experimental limits. 
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Figure 19: Predictions for the production cross sections times branching ratios of neutral elec- 
troweak vector resonances decaying to WW or dilepton final states in all three 
models. Only points with narrow resonances (P/m < 0.05) are shown. The dashed 
and solid curves show the expected and observed 95% C.L. experimental limits. 
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Figure 20: Predictions for the production cross sections times branching ratios of neutral elec- 
troweak vector resonances decaying to WW vs. dilepton final states in all three 
models. Only points with narrow resonances (P/m < 0.05) are shown. 

light quark pairs, with a branching ratio up to 40%, while the tt and bb hnal states have 
branching ratios below 10% each. The relative width of pc is around 10-50% when the 
decays to SM states are relevant, with closer to the upper end of this range. 

Charged resonances Among the three charged resonances, the lighter two are always 
nearly degenerate, with the lighter one being mostly the and the heavier one mostly 
the while the third charged resonance can be heavier and is mostly the axial vector 
resonance a^. The most important state for collider phenomenology is the second one 
since it is the only one with an appreciable Drell-Yan production cross section. Since its 
couplings to SM quarks are even weaker compared to the gluon resonance, it typically 
decays to quark partners, if kinematically allowed. If not, it decays to WZ and Wh 
with roughly equal branching ratios (cf. the discussion in sections 3.3.3 and 4.5.2). The 
branching ratio into tb is typically small but can reach 20% in corners of the parameter 
space. The branching ratio to iv is always below a percent. The other two states could 
be produced via vector boson fusion that we have neglected in our analysis since it is 
expected to be very small at the LHC (see [177] for a recent discussion). We note that 
the axial vector resonance typically decays to WZ and Wh with the largest branching 
ratios and we hnd BR(p 3 —>■ W^Z) ss 3BR(p^ —)• W^h). 

Neutral electroweak resonances Among the hve neutral uncoloured resonances, 
the two heaviest are usually mostly the axial vector resonances that have a small pro- 

^®Note that this means that in our numerical analysis, there are effectively no direct bounds on pc due 
to our requirement of a narrow width in the LHC searches, see section 3.3.2. 
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Figure 21: Predictions for the production cross sections times branching ratios of neutral elec- 
troweak vector resonances decaying to top quarks in all three models. Only points 
with narrow resonances (P/m < 0.05) are shown. The dashed and solid curves show 
the expected and observed 95% C.L. experimental limits. 


duction cross section in quark-antiquark collisions and (if produced via vector boson 
fusion) would decay with the largest branching ratios to WW and Zh with BR(p 4 5 —)> 
W^W^) ~ 3 BR (/94 5 —>• zh). Concerning the other three resonances, which are lin¬ 
ear combinations of the p^, p^, and they again preferably decay to a pair of 
quark partners. If this is kinematically disfavoured, they can decay to pairs of SM 
quarks, leptons, or W bosons, or to Zh. In the latter two cases, one typically has 
BR(/ 9 ° —^ W^W^) « BR(p? —^ Zh), as expected for an SU{2)l triplet (cf. the dis¬ 
cussion in sections 3.3.3 and 4.5.2). We find the branching ratios into electron or muon 
pairs to always be below 2 %, which can however be overcome by the higher experimental 
sensitivity, cf. figure 19 bottom. The branching ratio to light jets can be up to 30% in 
the U{2)^ models and up to 70% in U{3)^fj, while the one to bb can be up to 40% in 
all models. In the U{2)^ models, the decay to ti can come close to 100% in parts of the 
parameter space. 

In sections 4.3.5 and 4.4.4, we have already discussed that sizable NP contributions 
to the rare B decay Wilson coefficient Cg^, as is required if one wants to solve the 
anomalies in B physics discussed in section 3.2.5 in terms of new physics, requires a 
light, narrow neutral vector resonance with a large branching ratio to tt. In figure 21, 
we show the predictions for the production cross section times branching ratio of the 
neutral electroweak resonances decaying to ti at LHC with ^/s = 8 TeV, compared 
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to existing ATLAS and CMS analyses. As in the previous plots, we only show points 
with narrow resonances (T/m < 0.05). At masses below 1 TeV, the U{2)‘^q points 
correspond to the ones generating sizable NP effects in Cg^. The points in right-handed 
compositeness only start at higher masses because the relative width is typically larger 
than 5%. The plot shows that cross sections not far from what LHC has probed in run 
1 are attainable in all models. We conclude that this channel remains a promising probe 
at run 2, and discoveries are possible both for low and high masses. 

5. Summary 

In this paper, we have performed a comprehensive numerical analysis of a four-dimen¬ 
sional pNGB Higgs model based on the symmetry breaking coset SO{b)/SO{A) with 
quark partners transforming as fundamentals of 50(5). The model features a calculable 
one-loop Higgs potential and a custodial protection of the ZbihL coupling. We have 
included constraints from electroweak precision tests, flavour physics, Higgs production 
and decay, contact interaction searches, as well as direct searches for quark and vector 
resonances. We have considered three different flavour symmetries, all of them exact 
in the composite sector and broken only by the composite-elementary mixing terms, 
namely U{2)^ or U{3)^ with left- or right-handed compositeness. Below, we summarize 
our main findings. 

• Model-independently, we have pointed out that there are holes in existing experi¬ 
mental searches for quark partners decaying to W or Z plus a top or bottom quark, 
particularly for quark partners around 350 GeV which are not covered by Tevatron 
or LHC searches, see figure 2. We call on the experimental collaborations to close 
these holes by reanalyzing existing data. Quark partners decaying to a boson and 
a light quark are still weakly constrained. 

• In our numerical analysis, we have not found a single valid parameter point for the 
f7(3)Lc flavour structure. Although not a formal proof, we think this is a strong 
indication that this flavour structure is not compatible with electroweak precision 
tests and radiative EWSB in the model setup considered by us. 

• We have shown that the three other flavour structures can be made compatible 
with all relevant constraints with a fine tuning Abg ^ 100, see figure 3. 

• We have demonstrated that first-row CKM unitarity is the most sensitive probe 
of light-quark compositeness in U{2)f^Q, while in right-handed compositeness the 
dijet angular distribution is most sensitive to it, cf. figures 4 and 9. 

• Higgs signal strengths are the cleanest observables to constrain the pNGB decay 
constant /, with small corrections due to light-quark compositeness, cf. figures 5 
and 10. 

• In meson-antimeson mixing in the U{2)^ models, the relations between and 
Bg mixing that are expected from a leading-order spurion analysis are strongly 
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violated for some of the valid points by terms that are formally of higher order in 
the spurion expansion. 

• In both U{2)^ models, all observables in B, Bg and K mixing can saturate their 
current experimental limits, while in U{3)^, this is true for the mass differences 
and €k, while CP violation in the B and Bg systems is SM-like. The best means 
to experimentally distinguish the models based on AF = 2 observables alone can 
be read off figures 6, 11, and 12: 

— In U{2)f^Q, the relative NP effect in ex compared to AM^ is always smaller (or 
equal), in U{2)^f^ it is always larger (or equal), while in U{3)^fj it is always 
equal. 

— In U{2)f^(j, there can be large NP effects in (pg which are typically correlated 
with an equal effect in <pd', in U{2)^q NP effects can be only in (j)s and not in 
4)^; in U{3)^q both phases are free from NP. 

• In the U{2)^ models, CP violation in D^-D^ mixing is small compared to the cur¬ 
rent experimental sensitivity, but could become relevant if the sensitivity improves 
by an order of magnitude. 

• The FCNC top decay t —?• cZ can reach a branching ratio of up to 10“® in U{2)^q 
but is negligible in the other models. 

• Rare B decays of the type b —)• si^i~ can not only receive Z-mediated contribu¬ 
tions, but also resonance-mediated contributions that can affect the Wilson coeffi¬ 
cient Cg^ which is required for a NP explanation of various anomalies in B physics 
data. These anomalies can be explained in all three models. In [/(3 )^q, NP affects 
the Wilson coefficient Cg® but not C^q. 

• Explaining the B physics anomalies implies the presence of a narrow neutral vector 
resonance around 1 TeV with a sizable branching ratio into ti or dijets with a 
production cross section just below what has been excluded in LHC run 1. 

• Various excesses in diboson events at a mass of roughly 2 TeV observed by ATLAS 
and CMS can be explained in all three models as well by the decay of a 2 TeV 
vector resonance, see figures 18 and 19. The solution possibly, but not necessarily, 
predicts a signal in dilepton events around the same mass as well, see figure 20. 

While we have limited ourselves to a single model with four different flavour structures 
in this work, there are several ways how our analysis could be generalized, such as 
studying non-minimal cosets, non-minimal couplings, different fermion representations, 
or different flavour structures, including more radical changes like disposing of partial 
compositeness for the first two generation quarks [57, 58] . It would also be interesting 
to include a more realistic lepton sector. Finally, a more accurate treatment of the top 
quark mass, of loop corrections to the ZbihL coupling, of renormalization group effects 
on FCNC operators, and of LHC constraints on singly produced fermion resonances 
would be very interesting to further scrutinize composite Higgs models in the future. 
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A. SO(5) conventions 

For concreteness, we will present the conventions for SO(5) generators and embeddings 
used by us. 

The group 50(5) can locally be expressed as 50(5) = SU{2)lXSU{2)rXSO{ b)/SO{A). 
Therefore, its 10 generators can be grouped into “left”, “right” and “coset”; 

’ ( 78 ) 

(T“«),, = - 5]5i) - , (79) 

m 

Then, 50(4) singlets 5 and bidoublets (with SU{2)lX SU{3)r quantum numbers 

(ni,n 2 )) can be embedded into 50(5) fundamentals via 

f \ ( Q^^ + Q~~ \ 

.. iQ++ - iQ~~ 

= Q+--Q-+ . (81) 

ig+- + ig-+ 

V gv V v/25 / 

B. Mass matrices 

In this appendix we give the expressions for the mass mixing matrices that were obtained 
in the M4DCHM5. 

B.l. Boson sector 

The pNGB structure of the M4DCHM Lagrangian leads to mixings between the elemen¬ 
tary and composite vector bosons of equal charge. In particular, the composite triplets 
and as well as the axial resonances will have neutral and charged compo¬ 
nents mixing with the elementary and gauge bosons. In addition, the neutral 
components will also mix with the U{l)x resonance p^. 
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For the neutral and charged vector bosons one finds the following mass matrices given 
in tab. 4. 

By the explicit mixing introduced in the Lagrangian one finds the following mass 
matrices for the gluon and their composite resonances. By construction this does not 
spoil invariance under the SM SU{3)c, which survives as a linear combination of the 
elementary and composite SU{3) symmetries as can be seen from the fact that the 
gluon mass matrices exhibit a massless eigenvalue. 



(84) 


B.2. Fermion sector 

After EWSB the elementary quarks mix with all resonances carrying the same elec¬ 
tric charge. By using the embedding (81) we express the components of the bidou¬ 
blet resonances in such a way that they have definite quantum numbers under the 
50(4) = SU{2)l X SU{2)ii custodial symmetry. 

For the up- and down-type quarks we find the mass matrices given in tab. 5. 

For the exotically charged fermion resonances the mass matrices are independent of 
the Higgs held. Thus, they do not give a contribution to the Higgs potential, which is 
clear since they do not mix with elementary helds. 



Of course, the helds used above still carry havour indices. As a consequence of this, all 
the entries tin the fermionic mass matrices actually are 3x3 matrices in havour space, 
promoting the up- and down-type mass matrices to 27 x 27 objects. The explicit form 
of the entries is model dependent and will be given in appendix C. 

Since we took the leptons as purely elementary, their mass matrices are just diagonal 
taking the SM values. 

C. Explicit form of the composite-elementary mixings 

In this appendix, we give the explicit havour structure of the composite-elementary 
mixings in the havour symmetric models. We use bases where all unphysical parameters 
have already been rotated away and all phases have been made explicit. 

^°For example, the field Qt~ is part of a composite bidoublet resonance with = | and it has 
eigenvalues and — | under and respectively. 
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■^^Boson,neutral 



Table 4: Mass matrices for the neutral and singly charged bosons in the M4DCHM5 



Table 5: Mass matrices for the up- and down-type fermions in the M4DCHM5 


I 


cos 




In U{3)lc, 


Au,, = Au 1 


^dL = ^Lb 1 , 


^ ^Ru 

At ^ = yt ( Ar, 


^Rt/ 


At = 

dR 


^Ajid 


Ars 


^Rbj 


Here, V is the CKM matrix with 3 angles and 1 phase. 
In U{3)ic, 


Au^ = Ft ( Ai, 


At ^ = Ant 1 


Aity 


^^Ld 

^dL = I Ais 


A 


^dR - ^Rb 1 • 


Lbj 


(87) 

( 88 ) 


(89) 

(90) 


In Ui2)lc, 


Aur — 


A 


Lu I ) ^UR 

Alt 


'Aid 

AdR = I Aid 


A 


Lb) 


In t/(2)|c> 


CuAru 

AL = ( s„AK„e-*“ 


^dR ~ SdARde CdARs edARbe 


-s„ARce*"“ \ 

CuArc e^Aflte^'t’* , 

Am J 
(91) 

/ CdARd -SdARse'-^d 


^dR 


itj>h 


\ 


A 


Rb 


(92) 


^u^Lu Su^Lc^ 


^ut — I Su^Lu^ 


CuAlc euAite^'^^ , At^ = 

Au J 




^A 


Ru 


A 


Ru 


Arij 


(93) 


/ CdAid — 

AdR = SdALde~^°‘^ 

\ Au 


CdAu erfAifeC*'^'’ , At^ = 


''ARd 


A 


Rd 


ArL) 


(94) 
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D. Constraints from the dijet angular distribution 


As discussed in section 3.2.6, experimental analyses of contact interactions typically only 
quote constraints on a single operator - or for individual operators, but only allowing 
one at a time. To correctly treat the case with simultaneous contributions from multiple 
operators, we follow the procedure outlined in [119]. In this paper, analytical expressions 
are given for the dijet cross section in bins of the dijet mass rrijj and the rapidity x- 
The most recent ATLAS and CMS analyses use multivariate techniques rather than 
considering only a ratio of bins. In our numerical analysis, we have thus adopted the 
following procedure: 


1. We identify the most sensitive bin in the experimental analysis; 

2. We compute the numerical coefficients V and Q dehned as in [119] for the 8 TeV 
LHC in the respective bin. 

3. We compute the NP contribution of all operators to the cross section in this bin; 

4. We multiply our result by an overall factor to exactly reproduce the 95% C.L. 
bound on the Wilson coefficient c'qq quoted in the experimental paper. 


In this way, our approximation of computing the cross section analytically and pretend¬ 
ing that only a single bin is relevant is only used for the relative contributions of the 
individual operators, while any overall change (such as k-factors) cancels out since we 
normalize to the bound obtained for the c^qq coefficient by the experimentalists. We 
have checked that the relative contributions are not very sensitive to changes in the bin 
chosen in the first step. 

In our numerical analysis, we use the bound from the most recent ATLAS analysis 
[178]. We assume the most sensitive bin to be the one with x < 3.32, rrijj > 3.2 TeV. 
We can then write the new physics contribution to the dijet cross section in this bin, 


(j^. = 

3] 


, da{pp 

dx - 


n) 


<1 


dx 


NP 


mjq>Z.2 TeV 


(95) 


normalized to the 95% C.L. cross section on this quantity extracted by reproducing the 
bound on c^qq quoted in [178], as 


a 


30 


{^jj)95% C.L. 


= 

A2 




(96) 


where A and B are given in eq. (16) of [119] and V', Q' are equal up to normalization 
to "P, Q defined in [119]. Numerically, we find 

V' = (0.36P'„ 0.12P',, 0.36P^^, 0.12^^^, 0.17PU OAdP^^) , (97) 

Q' = (0.013gL, 0.0069Q;,, 0.013g:^^, 0.0069g:,^, 0.0024g;^, 0.00097g;^) , (98) 
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where 

Pl = (4.93TeV)2, = (1.46 TeV)^ , P', = (3.82 TeV)^ . (99) 

QL = (7.93 TeV)^ , = (4.28 TeV)^ , = (6.95 TeV)^ . (100) 
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